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ABSTRACT 

Optical observations of the Type la supernova (SN la) 2005bl in NGC 4070, obtained 
from —6 to +66 d with respect to the B-band maximum, are presented. The photomet- 
ric evolution is characterised by rapidly-declining light curves (Ami5(i3)truc — 1-93) 
and red colours at peak and soon thereafter. With — —17.24 the SN is an 

underluminous SN la, similar to the peculiar SNe 1991bg and 1999by. This similar- 
ity also holds for the spectroscopic appearance, the only remarkable difference being 
the likely presence of carbon in pre-maximum spectra of SN 2005bl. A comparison 
study among underluminous SNe la is performed, based on a number of spectropho- 
tometric parameters. Previously reported correlations of the light-curve decline rate 
with peak luminosity and 7?,(Si) are confirmed, and a large range of post-maximum 
Sin A6355 velocity gradients is encountered. ID synthetic spectra for SN 2005bl are 
presented, which confirm the presence of carbon and suggest an overall low burning 
efficiency with a significant amount of leftover unburned material. Also, the Fe con- 
tent in pre-maximum spectra is very low, which may point to a low metallicity of the 
precursor. Implications for possible progenitor scenarios of underluminous SNe la are 
briefly discussed. 

Key words: supernovae: general - supernovae: individual: SN 2005bl - supernovae: 
individual: SN 1991bg - supernovae: individual: SN 1999by - supernovae: individual: 
SN 1998de - galaxies: individual: NGC 4070. 



1 INTRODUCTION 



* Based on observations at ESO-Paranal, Prog. 075.D-0662(B) The history of underluminous SNe la is a typical example of 

t E-mail: tauben@mpa-garching.mpg.de the ever-recurring pattern in which knowledge about nature 
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is accumulated. Usually, in the beginning there is the obser- 
vation of a phenomenon, followed by a successful theoretical 
explanation. However, as further experiments or observa- 
tions are carried out in order to confirm the newly developed 
theoretical ideas, often an ever higher degree of diversity and 
ever more exceptions from the simple rules are found the 
closer the subject of interest is studied. The need for refined 
and more complex theories to obtain a realistic description 
of the involved processes becomes evident. 

In the case of SNe la, first a class of cosmic ex- 
plosions apparently similar in absolute luminosity ("stan- 
dard candles") and spectroscopic appearance was iden- 
tified. These events were explained as the disruptions 
of white dwarfs which had accreted matter until they 
reached their stability limit close to the Chandrasekhar 
mass (Afch). However, in 1991 the paradigm of SN la 
homogeneity had to be relaxed a lot. This was triggered 
by the observation of two peculiar SNe la, which there- 
after served as prototypes of newly-defined SN la subclasses 
with distinct spectrophotome t ric pro p erties. One of these , 
SN 1991T (iFilippenko et all [l992al: iPhillips et al.l Il992l: 



Ruiz-Lapuente et al.l 19921 : iMazzali. Danziger fc Turattd 
19951 ). was UP to 0.6 mag brighter than average SNe la. 



and characterised by a hot early-time spectrum with 
strong Felll features and weak or absent Si ll and 
Sii lines. The other one, SN 1991bg (IFilippenko et al.l 
1992bl: iLeibundgut et all Il993l: iRuiz-Lapuente et all 1199^ 



Turatto et al.lll996l : iMazzali et al. 1 19971 ) was even more de- 
viant, with low ejecta velocities and a cool spectrum domi- 
nated by intermediate-mass-element (IME) lines and partic- 
ularly strong Ol and Till. Moreover, it had unusually red 
colours at early phases, and was underluminous by about 
2 mag at peak (hereafter we will refer to such appearance 
as 91bg-like). Hence, quasi instantaneously the luminosity 
range of SNe la had increased to a factor of ten between the 
brightest and the faintest objects, proving th at they were no 
standard candles. However, two years later IPhillipsI (Il993l ) 
realised a tight correlation between peak luminosity and de- 
cline rate in the B band. Thi s relation and revised versions 
of it (e.g. iPhiUips et al.lll999l ) turned SNe la into standard- 
isable candles, and hence made them an extremely useful 
tool for precision cosmologjQ. 

In the following years, several more 91bg-like SNe la 
were discovered, but the available data set grew much less 
rapidly than for ordinary SNe la. From the results of the 
Lick Observatory Supernova Search (LOSS) and the Bei- 
jing Astronom ical Observatory Supernova Survey (BAOSS), 
iLi et all (|200ll ) estimated that about 16 % of all SNe la are 
of the 91bg-like variety. This may still be an underestimate, 
as their low intrinsic luminosity makes 91 b g-like SNe prone 
to Malmquist bias; nevertheless iLi et al.l l|200ll ') estimated 
this effect to be negligible in thei r sample. Statistical stud- 
ies l|Hamuv et al.lll996al . [2000l :|Howell 2001) have shown that 
SNe la occur in all host-galaxy types, but revealed a corre- 
lation between SN decline rate and host morphology, with 



a clear tendency for 91bg-like SNe to be associated with 
early-type hosts and hence old stellar populations. 

While the single-degenerate (SD) Chandrasekhar-mass 
model has survived as the favoured scenario for the normal 
and indeed rather homogeneous SNe la, a number of alter- 
native models have been suggested for the 91bg-like sub- 
class. Ideas include double-degenerate (DD) explosions of 
merging white dwarfs, sub-Chandrasekhar-mass explosions 
triggered by detonation of th e accreted helium layer (cf. 
iHillebrandt fc Niemeveij |2000| . for a review), and defiagra- 
tions in strongly rotating white dwarfs, where the turbulent 
propagation of the flame fro nt is suppressed by the differen- 
tial rotation (|Pfannesll2006l ). Still, the notion that 91bg-like 
SNe are - in terms of the underlying explosion model - no 
different from ordinary SNe la, and that the only discrimi- 
nating parameter is the mass of synthesised ^®Ni, has sup- 
porters in the SN la community. No conclusive evidence for 
any of these ideas has been found so far. 

In this paper we present the joint data set of SN 2005bl 
obtained by the European Supernova Collaboration (ESC jf| 
and the Carnegie Supernova Project (CSPfl Since these ob- 
servations are among the earliest ever obtained for a 91bg- 
like SN, they may help to better constrain possible progeni- 
tor and explosion models. The observations and techniques 
applied for data reduction and calibration are discussed in 
Section[2l In Section[3]we estimate the distance of SN 2005bl 
and the extinction along the line of sight. Sections |4] and [5] 
are devoted to the analysis of the light curves and spectra, 
respectively. Results of ID spectrum synthesis calculations 
are presented in Section [B] and a comparison with other un- 
derluminous SNe la is performed in Section [T] where we also 
discuss the impact of SN 2005bl on our picture of SN la ex- 
plosions. A short summary of the main results is given in 
Section H 



2 OBSERVATIONS AND DATA REDUCTION 

SN 2005bl (z = 0.024) was discovered in the course of 
the Lick Observatory Supernova Search programme (LOSS) 
with the Katzman Automatic Imaging Telescope (KAIT) on 
UT 2005 April 14.34 and 15.36 at unfiltered magnitudes of 
18.8 and 18.3, respectively (Shimasaki & Li 2005). The SN 
was not detected on images obtained with the same setup on 
UT 2005 March 11.33 to a limiting magnitude of 19.5. Based 
on spectra taken with the Las Campanas 2.5 m du Pont 
Telescope (-1- WFCCD spectrograph) and the Fred Lawrence 
Whipple Observatory 1.5 m Telescope (-I-FAST), SN 2005bl 
was classified as SN la, probably belonging to the 91bg- 
like variety given the similarity of the spectra with those of 
SN 1999by a few days before maximum light l|Morrell et al.l 
12005': 'Gallagh er et al.|[2"005ah . The SN is located in the ellip- 
tical galaxy NGC 4070 (de Vaucouleurs morphological type 
-4.9; LEDAQ), projected on a region of steep yet smooth 
background variation (Fig. [TJ . Details on the SN and host 
galaxy properties are summarised in Table [1] 



^ Several years later ev en more peculiar objects such as 
SN 2002CX llLi et al.ll2003t) were discovered which did not obey 
the Phillips relation, showing that not all SNe la are standardis- 
able. 



^ http: // www. mpa-garching. mpg.de /~rtn/ 
^ http://www.cspl.lco.cl/~cspuserl/PUB/CSP.html 
^ Lyon-Meudon Extragalactic Database 
http: //leda. univ-lyonl.fr/ 
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Figure 1. /?-band image of the SN 2005bl field taken with the 
Calar Alto 2.2 m Telescope + CAFOS on UT 2005 May 14. The 
field of view is 9 X 9 arcmin'^, and the local sequence stars are 
indicated. In the upper-left corner a twice-enlarged blow-up of 
the SN and its host galaxy is shown. 

Although SN 2005bl was too distant to fulfil the formal 
selection criteria of the ESC, optical follow-up observations 
were performed owing to the peculiarities found in our first 
spectrum, obtained almost at the same time as the classifi- 
cation spectra by other groups. Ifowever, given the dimness 
of the SN, an intensive coverage as for other ESC targets was 
out of reach. In particular, the early part of the light curves 
was not well sampled owing to bad weather and scheduling 
constraints. At the same time, the CSP collaboration started 
their follow-up of SN 2005bl, focussing mainly on the photo- 
metric evolution near maximum light. Hence, the two data 
sets were almost perfectly complementary. 

2.1 Photometric data 

Optical photometry of SN 2005bl was acquired from one 
week before to about two months after maximum light in B. 
The basic data reduction (bias subtraction, overscan correc- 
tion and flat-fielding) was perfo r med using sta ndard routines 
in lRAlPI (|Massev fc Davislll992l : lMassevlll997l '). The local se- 
quence of stars in the SN fie ld shown in Fig . [T] was calibrated 
with r espect to a number of lLandoltl (|l992l . for UBVRI) and 
Sloan (jSmith et al.ll2002l . for ugriz) standard fields on sev- 
eral photometric nights. The magnitudes of the calibrated 
local sequence, listed in Tabled were used subsequently to 
determine the SN magnitudes in a relative measurement. 

5 NASA/IPAC Extragalactic Database 

http: //nedwww. ipac.caltech.edu/ 
^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc, under contract to the National 
Science Foundation. 



Table 1. Properties of SN 2005bl and its host galaxy. 



NGC 4059 / NGC 4070° 



a 


12h04mll?43 


1 


5 


-1-20° 24' 37'.' 7 


1 


redshift 


0.02406 ± 0.00008 


2 


recession velocity v 


7213 ± 24 kms-l 


2 


Virgo 


7330 ± 28 kms-i 


2 


"CMB 


7534 ±33 kms-i 


2 


distance modulus /x'' 


35.10 ±0.09 mag 


2 


apparent corr. B magnitude 


13.75 ±0.15 


1 


morphological type*^ 


E, -4.9 


1 


Galactic reddening E{B — V) 


0.028 mag 


3 


SN 2005bl 


a 


12h04ml2!32 


4 


5 


+20°24'24'.'8 


4 


offset from galaxy centre 


13'.'9 E, 11'.'2 S 


4 


host reddening E(B — V) 


0.17 ±0.08 mag 


5 


Ami5(B)truo 


1.93 ±0.10 


5 


JDmax in U 


2 453 481. 4 ±0.3 


5 


JDmax in B 


2 453 482. 6 ±0.3 


5 


JDmax in V 


2 453 484.9 ±0.3 


5 


JDmax in R 


2 453 485. 9 ±0.3 


5 


JDmax in / 


2 453 487.0 ±0.3 


5 


JDmax in g 


2 453 483.0 ±0.3 


5 


JDmax in 2 


2 453 487.0 ±3.0 


5 


^^max 


19.14 ±0.18 


5 


^max 


18.68 ±0.04 


5 


Vmax 


17.87 ±0.03 


5 


-Rmax 


17.55 ±0.03 


5 


-^max 


17.38 ±0.04 


5 


5max 


18.20 ±0.07 


5 


2max 


17.77 ±0.12 


5 


M(7_max 


-16.91 ±0.43 


5 


Mb, max 


-17.24 ±0.34 


5 


My^max 


-17.85 ±0.27 


5 


Mjj^max 


-18.06 ±0.23 


5 


M/,max 


-18.10 ±0.18 


5 


Mj^^max 


-17.55 ±0.29 


5 


,max 


-17.67 ±0.20 


5 



" The galaxy is listed twice in the NGC catalogue. 
V corrected for Local-Group infall onto Virgo cluster 
V corrected to the CMB reference frame 
from VQMB, using Hq = 72kms~^Mpc~^ 
numerical code according to de Vaucouleurs 

1: LEDA; 2: NElfl; 3:[ Schlegcl. Finkbe iner fc Da^ll998l : 4: 
jPuckett fc Laneoussisii2005i : 5: this work 



For the sequence-star magnitudes in the Sloan ugri bands, 
the reader is referred to Contreras et al. (in prep.). 

Although the host-galaxy background at the projected 
SN site seemed to be smooth , wc ap p lied t he template- 
subtraction technique (jFilippenko et aU Il986l ) in order to 
eliminate any possible contamination from the host galaxy, 
a justified concern given the faintness of the SN. For this pur- 
pose we acquired templates in BVRI with CAFOS mounted 
on the Calar Aho 2.2 m Telescope on UT 2006 March 28, and 
in uBVgri with the Las Campanas 2.5 m du Pont Telescope 
from UT 2006 April 3 to 30, about one year after the ex- 
plosion when the SN had faded from visibility. For the ESC 
data, the galaxy subtraction was performed using the iraf 
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Table 2. Magnitudes of the local sequence stars in the field of SN 2005bl (Fig. [TJ. 



ID 



B 



V 



R 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 



15.207 ± 0.022 
16.387 ± 0.022 
18.340 ± 0.019 



19.274 
19.538 
18.656 



0.017 
0.022 
0.049 



18.355 ± 0.049 



14.786 
15.931 
18.633 
18.154 
17.464 
19.132 
19.488 
19.396 
18.924 
20.035 
18.267 
18.346 
18.278 
19.083 



± 0.013 
± 0.012 
± 0.021 
± 0.014 
± 0.009 
± 0.014 
± 0.020 
± 0.021 
± 0.016 
± 0.020 
± 0.013 
± 0.014 
± 0.011 
± 0.013 



13.952 
15.048 
18.110 
16.658 
16.962 
18.220 
18.979 
18.659 
18.454 
18.640 
17.539 
17.873 
17.206 
18.198 



± 0.014 
± 0.014 
lb 0.016 
± 0.020 
± 0.027 
± 0.015 
± 0.014 
± 0.019 
± 0.018 
± 0.022 
± 0.021 
± 0.020 
± 0.017 
± 0.023 



13.502 
14.510 
17.767 
15.665 
16.621 
17.690 
18.629 
18.251 
18.115 
17.532 
17.104 
17.510 
16.509 
17.666 



± 0.010 
± 0.011 
± 0.013 
± 0.021 
± 0.007 
± 0.018 
± 0.008 
± 0.013 
± 0.020 
± 0.014 
± 0.015 
± 0.016 
± 0.018 
± 0.013 



13.080 
14.011 
17.385 
14.792 
16.269 
17.247 
18.272 
17.824 
17.834 
16.260 
16.687 
17.196 
15.893 
17.193 



0.013 
0.012 
0.012 
0.017 
0.012 
0.010 



± 0.015 



0.017 
0.018 
0.022 
0.011 
0.011 
0.025 
0.020 



13.386 ± 0.013 
14.295 ± 0.024 
17.710 ± 0.074 



18.526 ± 0.035 
18.096 ± 0.030 
18.112 ± 0.040 

16.965 ± 0.014 



plug-in SVSUB written by S.V. (based on isis). The instru- 
mental SN magnitudes were determined in the background- 
subtracted images with point-spread function (PSF) fitting 
photometry using the software package SNOOPY, specifically 
designed for this purpose by F. Patat and implemented in 
IRAF by E. Cappellaro. For the z band no templates were 
available, and the measurements were performed with or- 
dinary background-fitting PSF photometry in SNOOPY. A 
check in the BVRI bands showed a good agreement between 
the two methods. Therefore, we are confident that also our 
a-band photometry is sufficiently reliable. 

The calibration of the SN magnitudes to the de- 
sired standard photometric systems, the Johnson/Cousins 
system (J Bcsscll 1990) for UBVRI and the Sloan system 
l|Fukugita et al.lll996l ) for g and z, was complicated by the 
variety of filters mounted at the various telescopes, some of 
which deviated strongly from the standard. Fig.[2]shows the 
BVRI and ri response curves of all instruments used for the 
SN follow-up. 



Bessell photometry 

In order to compensate for the differences in the trans- 
mission curves, and to report the magnitudes on the 
Bessell system, we made use of the "^'-correction " tech - 
nique based o n the p rescription of IStritzinger et"al] l|2002h . 
iPignata et al.l (|200J), and references therein. S-corrections 
were computed for the BVRI and Sloan ri data, based on 
our spectra of SN 2005bl. Unfortunately, the spectra do 
not fully cover the [/-band region, so that no ?7-band S- 
corr ection could be derived. Instead, transformation formu - 
lae (IJordi. Grebel fc Ammonll2006l : IZhao fc Newbergj|2006l ) 
were employed to convert the Swope M-band data to the 
Bessell system. Similarly, no S-correction was applied to 
data later than -1-50 d, since our spectroscopic follow-up ends 
already well before. The strong deviation of the MONICA 
V filter (actually a Roeser BV filter) from the Bessell de- 
scription (see Fig. [2} , along with some lacking information 
required to reconstruct the response curve, makes the MON- 
ICA l/-band correction less reliable. 

The redshift of SN 2005bl is not negfigible. Therefore, 
in computing the S-correction the spectra were shifted in 
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4000 5000 6000 7000 8000 9000 10000 
wavelength [A] 

Figure 2. Instrumental BVRI / ri passbands use d for the obser- 
vations of SN 2005bl. The standard Sesselll l|l990h curves are also 
displayed in the figure for comparison. 



order to account also for the if-correction. Since the same 
restrictions as before apply also here, no isT-correction was 
performed for the U band and for any data later than -1-50 d. 

In Table[3]the fully calibrated and - whenever possible - 
S- and Tf-corrected UBVRI Bessell magnitudes of SN 2005bl 
are reported, together with their uncertainties. Both the in- 
trinsic Bessell data and the transformed Sloan photometry 
entered into this Table. Table |3] shows the combined 5*- and 
Tf-correction, i.e., the difference between the magnitudes of 
Table |3] and those obtained with a first-order colour-term 
calibration. The differences are significant, so that the lack 
of the S- and isT-correction for U may introduce considerable 
uncertainties. 
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Table 3. S- and /S'-corrected Bessell magnitudes of SN 2005bl.° 



JD^ 


Epoch*^ 


U 


B 








V 






R 








/ 




Telescope 


Seeing'* 


476.61 


—6.0 


1 Q c;qi _|_ n 1 1 4 


19.364 


± 


0. 


,022 


18.950 


-|- 


0.030 


18.776 


-|- 


0. 


,024 


18.599 


-|- 


0.045 


SWO 


1.41 


477.63 


—5.0 


1 Q "^^^Q -1-0117 


19.135 


± 


0. 


,022 


18.695 


-|- 


0.029 


18.523 


-|- 


0. 


,023 


18.402 


-|- 


0.038 


SWO 


1.34 


478.66 


—3.9 




18.955 


± 


0. 


022 


18.461 


-|- 


0.029 


18 290 


-|- 


0. 


022 


18.120 


-|- 


0.036 


SWO 


1.68 


479.64 


—3.0 




18.818 


± 


0. 


,029 


18.307 


-|- 


0.031 


18.116 


-|- 


0. 


,030 


17.948 


-|- 


0.043 


SWO 


1.57 


479.68 


—2.9 












18.393 


-|- 


0.042 


18.120 


-|- 


Q 


,030 


17.883 


-|- 


0.034 


TNG 


1.47 


480.65 


— 1.9 


1 Q 1 r:i -1- 1 zl4 


18.766 


± 


Q 


,025 


18.150 


-|- 


0.029 


17.883 


-|- 


Q 


,028 


17.770 


-|- 


0.038 


SWO 


1.44 


481.61 


— 1.0 


1Q IOC 1 r\ 174 


18.709 


± 


Q 


,027 


18.052 


-|- 


0.030 


17.798 


-|- 


Q 


,022 


17.650 


-|- 


0.036 


SWO 


1.45 


483.62 


1.0 




18.682 


± 


Q 


,040 


17.899 


-|- 


0.030 


17.636 


-|- 


Q 


,026 


17.513 


-|- 


0.039 


SWO 


1.54 


484.65 


2.1 


1 Q 9Qfi -1- 9Q1 


18.733 


± 


Q 


,042 


17.839 


-|- 


0.031 


17.526 


-|- 


Q 


,043 


17.414 


-|- 


0.045 


SWO 


1.52 


485.39 


2.8 




18.872 


± 


Q 


,036 


17.887 


-|- 


0.051 


17.563 


-|- 


Q 


,035 


17.361 


-|- 


0.051 


LT 


0.84 


489.32 


6.7 












18.172 


-|- 


0.280 
















WD 


1 89 


489.59 


7.0 


90 94S -1- n 1 '^9 


19.472 


± 


Q 


031 


18.146 


-|- 


0.027 


17 692 


-|- 


Q 


021 


17.416 


-|- 


0.035 


SWO 


1.84 


490.52 


7.9 


90 fi4n -1- 1 '^'^ 


19.653 


± 


Q 


,026 


18.194 


-|- 


0.025 


17.781 


-|- 


Q 


,025 


17.435 


-|- 


0.044 


LT 


0.81 


4Qn fin 


8.0 


90 fi4n -1- 1 fi'^ 


19.687 


± 


Q 


,023 


18.258 


-|- 


0.028 


17.791 


-|- 


Q 


,022 


17.456 


-|- 


0.038 


SWO 


1.34 


491.51 


8.9 




19.824 


± 


Q 


,046 


18.296 


-|- 


0.048 


17.845 


-|- 


Q 


,027 








WD 


1.60 


491.57 


9.0 


90 fi9n -1- 1 40 


19.834 


± 


Q 


,022 


18.359 


-|- 


0.029 


17.836 


-|- 


Q 


,026 


17.477 


-|- 


0.039 


SWO 


1.37 


493.35 


10.8 












18.553 


-|- 


0.028 


17.959 


-|- 


Q 


,025 


17.503 


-|- 


0.031 


LOl 


2.41 


498.38 


15.8 


















18.551 


-|- 


Q 


,045 


17.805 


-|- 


0.052 


LT 


0.95 


500.51 


17.9 












19.258 


-|- 


0.184 


18.722 


-|- 


Q 


,041 








WD 


2.38 


ou^ .oo 


19.8 


















18.807 


-|- 


Q 


034 


18.192 


-|- 


0.063 


LT 


0.70 




20.8 




20.912 


± 





,073 


19.334 


-|- 


0.040 


18.960 


-|- 


0, 


,023 








WD 


1.55 


ouo. ou 


20.8 




20.893 


± 


Q 


,133 


19.429 


-|- 


0.046 


18.965 


-|- 


Q 


,074 


18.361 


-|- 


0.059 


Ekar 


2.57 


504.49 


21.9 


















19.035 


-|- 


0, 


,035 


18.325 


-|- 


0.057 


LT 


0.67 


ouo.4y 


oo o 
zz .y 




20.936 


± 


U 


,UOU 


ly .DUy 


It 


U.uoy 


101 fin 
ly . ioU 


_(_ 
It 


U. 


.U4o 


^ Q A lA 


It 


n c\c\A 

U.UD4 


i^ana 


i.zy 


506.41 


23.8 




21.052 


± 


0, 


,081 


19.597 


± 


0.037 


19.216 


it 


0, 


,028 


18.571 


± 


0.035 


Caha 


1.91 


510.53 


27.9 












19.693 


± 


0.119 
















WD 


1.70 


511.36 


28.8 












19.687 


± 


0.092 


19.530 


± 


0, 


,045 








WD 


1.70 


512.55 


29.9 




21.168 


± 


0, 


,176 


19.847 


± 


0.063 


19.480 


± 


0, 


.036 


18.874 




0.051 


SWO 


1.56 


516.48 


33.9 


















19.877 


± 


0, 


,073 








WD 


1.78 


518.43 


35.8 












20.044 


± 


0.068 


20.020 


± 


0, 


,042 








WD 


1.52 


519.40 


36.8 




21.421 


± 


0, 


.075 


20.031 


± 


0.042 


19.925 


± 


0, 


.055 


19.310 




0.096 


LT 


1.01 


520.44 


37.8 




21.407 


± 


0, 


.124 


20.205 


± 


0.055 


20.002 


± 


0, 


,045 








WD 


2.07 


521.39 


38.8 












20.122 


± 


0.040 


20.009 


± 


0, 


.036 


19.358 




0.094 


LT 


0.96 


524.42 


41.8 




21.379 


± 


0, 


,097 


20.141 


± 


0.069 


20.159 


± 


0, 


,071 








WD 


1.40 


538.38 


55.8 












20.685 


± 


0.168 
















WD 


1.75 


548.41 


65.8 












20.984 


± 


0.280 


21.183 


± 


0, 


,347 








WD 


2.15 



" No S- and /ST-correction applied in the f/ band and to any data after +50 d. ^ JD — 2 453 000.00 Epoch in days with respect to 
the B-band maximum JD 2 453 482.6 ±0.5. Average seeing in arcsec over all filter bands. 

SWO = Las Campanas 1.0m Swope Telescope + CCD; http;//virww.lco.cl/telescopes-information/henrietta-swope/ 
TNG = 3.58 m Telescopio Nazionale Galileo + DOLORES; http://www.tng.iac.es/instruments/lrs/ 
LT = 2.0 m Liverpool Telescope + RATCAM; http://telescope.livjm.ac.uk/Info/TelInst/Inst/RATCam/ 
WD = 0.8m Wendelstein Telescope + MONICA; http://www.wendelstein-observatorium.de/monica/monica_en.html 
LOI = 1.52 m Loiano Telescope + BFOSC; http://www.bo.astro.it/loiano/152cm.html 
Ekar = Asiago 1.82 m Telescope + AFOSC; http://www.oapd.inaf.it/asiago/2000/2300/2310.html 
Caha = Calar Alto 2.2m Telescope + CAFOS SiTe; http://www.caha.es/CAHA/Instruments/CAFOS/ 



Sloan photometry 



2.2 Spectroscopic data 



Besides transforming the Swope and Liverpool Sloan-filter 
observations to the Bessell system via S'-corrections, light 
curves in the Sloan photometric system itself were con- 
structed. The magnitudes were calibrated through first- 
order colour equations without employing S-corrections, 
since at both telescopes the filters are close to the standard 
Sloan prescription. Table [S] reports the Sloan griz photom- 
etry from the Liverpool telescope. For the original Swope 
photometry see Contreras et al. (in prep.). 



Details of the spectroscopic observations of SN 2005bl 
are reported in Table |6l All two-dimensional spectroscopic 
frames were first debiased and flat-fielded, before an op- 
timal , variance-weighted extraction of the spectra (|Hornel 
1 19861 : iMassev. Valdes fc Barnes! [l99^ ) was performed using 
the IRAF routine APALL. Wavelength calibration was accom- 
plished with the help of arc-lamp exposures or, whenever 
this was not possible, using the night-sky lines. The in- 
strumental response functions required for flux calibration 
were determined from observations of the spectrophotomet- 
ric standard stars reported in Table |6] Whenever no stan- 
dard had been observed, the sensitivity curve obtained on 



6 Taubenberger et al. 



Table 5. Sloan photometry of SN 2005bl obtained with the Liverpool telescope. No S- and X-correction has been applied. For the Swope 
data see Contreras et al. (in prep.). 



JD" 


Epoch'' 


g r 


i 


z 


485.4 


2.8 


18.295 ± 0.115 17.669 ± 0.023 


17.866 ± 0.051 


17.769 ± 0.080 


490.5 


7.9 


17.877 ± 0.025 


17.870 ± 0.031 


17.808 ± 0.039 


498.4 


15.8 


18.676 ± 0.046 


18.308 ± 0.092 


18.113 ± 0.080 


502.4 


19.8 


18.921 ± 0.039 


18.760 ± 0.059 




504.5 


21.9 


19.138 ± 0.047 


18.878 ± 0.077 


18.548 ± 0.095 


519.4 


36.8 


19.998 ± 0.071 


19.982 ± 0.093 


19.548 ± 0.165 


521.4 


38.8 


20.073 ± 0.043 


20.081 ± 0.100 


19.680 ± 0.240 



" JD - 2 453 000.0 Epoch in days with respect to B-band maximum (JD = 2 453 482.6 ± 0.5). 
Table 6. Journal of spectroscopic observations of SN 2005bl. 



UT Date JD" Epoch Exposure Airmass Tel. Grism Range [A] Res. [A]'= Standards 



05/04/16 


476.6 


-6.0 


900 s 


X 


3 


1.54 


DUP 


blue 


3800 


- 9200 


6 


L745-46A, LTT7987 


05/04/17 


477.6 


-5.0 


2400 s 






1.54 


Caha 


b200 


3500 


- 8800 


10 


BD-f 33 2642 


05/04/19 


479.6 


-3.0 


900 s 


X 


3 


1.54 


DUP 


blue 


3700 


- 9000 


6 


L745-46A, LTT7379 


05/04/19 


479.7 


-2.9 


1200 s 






2.06 


TNG 


LR-R 


5000 


- 9750 


11 


Hz44 


05/04/26 


487.4 


4.8 


1500 s 






1.09 


TNG 


LR-B 


3300 


- 8000 


13 


Feige 34 


05/04/26 


487.4 


4.8 


1500 s 






1.15 


TNG 


LR-R 


5000 


- 9750 


12 


Feige 34 


05/05/04 


495.5 


12.9 


1500 s 






1.04 


TNG 


LR-B 


3500 


- 8900 


14 


Feige 56 


05/05/04 


495.5 


12.9 


1500 s 






1.08 


TNG 


LR-R 


5000 


- 9750 


12 


Feige 56 


05/05/11 


502.4 


19.8 


1500 s 


X 


2 


1.01 


TNG 


LR-B 


3300 


- 8000 


14 


Feige 66 


05/05/11 


502.4 


19.8 


2700 s 


X 


2 


1.14 


Caha 


r200 


6200 


- 9750 


11 


Feige 34 


05/05/14 


505.4 


22.8 


2400 s 


X 


2 


1.11 


Caha 


b200 


3500 


- 8800 


14 


BD+33 2642 


05/05/24 


515.6 


33.0 


1800 s 






1.93 


VLT 


300V + GG375 


4200 


- 9600 


9 


LTT7987 



" JD - 2453 000.0 Relative to B-band maximum (JD = 2453482.6 ± 0.5). " Full-width at half maximum (FWHM) of isolated, 

unblended night-sky lines. 

DUP = Las Campanas 2.5 m du Pont Telescope + WFCCD; http://www.lco.cl/telescopes-information/irenee-du-pont/ 
Caha = Calar Alto 2.2m Telescope + CAFOS SiTe; http://www.caha.es/CAHA/Instruments/CAFOS/ 
TNG = 3.58m Telescopio Nazionale Galileo + DOLORES; http://www.tng.iac.es/instruments/lrs/ 
VLT = ESO 8.2 m Very Large Telescope UTl -|- F0RS2; http://www.cso.org/instruments/fors2/ 



a different night with the same instrumental configuration 
was used. Atmospheric extinction correction was applied us- 
ing tabulated extinction coefficients for each telescope site. 
Telluric features were identified in the spectra of the spec- 
trophotometric standard stars and removed from the SN 
spectra. To check the calibration, the spectroscopic fluxes 
were transformed in to magnitudes by integrating the spec- 
tra convolved with iBesselll (|l990l ) fllter functions. When- 
ever necessary, the spectral fluxes were adjusted to match 
the contemporaneous photometry. Finally, spectra of simi- 
lar quality obtained during the same night were combined 
to increase the signal-to-noise ratio (S/N); if the wavelength 
range of these spectra was different, they were averaged in 
their overlap region. 

As a consequence of charge-transfer-efficiency problems 
of the DOLORES CCD, it was not possible to remove the 
night-sky emission in the TNG spectra cleanly. A pattern of 
negative and positive residuals was left, sometimes at wave- 
lengths coinciding with spectral features of the SN, thus lim- 
iting the reliability of line-depth measurements. To mitigate 
this problem, the lines affected were fitted with polynomials, 
excluding from the fit the regions of strongest residuals, and 
then the depth was determined in the polynomial curves. 



DISTANCE, EXTINCTION AND HOST 
GALAXY PROPERTIES 



Like most 91bg-like SNe la (|GaUagher et al.1 l2005bh . 
SN 2005bl exploded in an early-type host, the elliptical 
galaxy NGC 4070. Other prominent examp les are SN 1991bg 
itself a nd SN 1997cn in ellipt i cal hosts (iFilippenko et aU 
Il992bl : iLeibundgut et al.l 1 19931 : iTuratto et al.Ml996l , llOoi T 



and SN 1998de w hose host galaxy was of type SO 
jModiaz et al.l l200l') PI Since early- type galaxies are assumed 
to have experienced no significant star formation over long 
times, this behaviour might be indicative of 91bg-like SNe la 
originating from an old stellar population. Even if one al- 
lows for some more recent star formation in elliptical galax- 
ies (e.g. triggered by mergers or interaction), the relative 
paucity of 91bg-like SNe in late-type, actively star-forming 
galaxies clearly disfavour young or intermediate-aged stellar 
progenitors or progenitor systems. 

A widely accepted paradigm for elliptical galaxies is 
that they have only little interstellar gas and dust (but see 



^ SN 1999by dCarnavich et al.ll2004l) was hosted by NGC 2841, 
a spiral galaxy, but with the s pectral appearance of an elliptical 
galaxy l lGallagher et al.ll2005bll . 
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Table 4. Amount of S- and X-correction contained in the mag- 
nitudes reported in Table |3] No 5- and /S'-correction has been 
applied to the (7-band data. 



JD" 


Epoch' 


B 


V 


R 


/ 


Tel.': 


476.6 


—6.0 


—0.020 


0.030 


0.033 


0.147 


swo 


477.6 


—5.0 


—0.038 


0.029 


0.044 


0.177 


swo 


478.7 


—3.9 


—0.028 


0.011 


0.045 


0.162 


swo 


479.6 


—3.0 


—0.054 


0.011 


0.044 


0.151 


swo 


479.7 


—2.9 




0.123 


0.193 


0.252 


TNG 


480.7 


-1.9 


-0.032 


0.007 


0.022 


0.093 


SWO 


481.6 


— 1.0 


—0.062 


—0.003 


0.043 


0.110 


SWO 


483.6 


1.0 


—0.078 


—0.044 


0.088 


0.102 


SWO 


484.7 


2.1 


—0.099 


—0.042 


0.048 


0.085 


SWO 


485.4 


2.8 


—0.077 


—0.002 


0.136 


—0.054 


LT 


489.3 


6.7 




0.137 






WD 


489.6 


7.0 


—0.113 


—0.039 


0.187 


0.060 


SWO 


490.5 


7.9 


—0.116 


—0.034 


0.166 


—0.042 


LT 


490.6 


8.0 


—0.128 


—0.065 


0.213 


0.049 


SWO 


491.5 


8.9 


—0.233 


—0.034 


0.074 




WD 


491.6 


9.0 


—0.142 


—0.074 


0.190 


0.035 


SWO 


493.4 


10.8 




—0.029 


0.027 


—0.018 


LOI 


498.4 


15.8 






0.126 


—0.105 


LT 


500.5 


17.9 




—0.064 


0.096 




WD 


502.4 


19.8 






0.096 


-0.209 


LT 


503.4 


20.8 


-0.237 


-0.065 


0.074 




WD 


503.4 


20.8 


-0.166 


-0.067 


0.055 


-0.017 


Ekar 


504.5 


21.9 






0.113 


-0.220 


LT 


505.5 


22.9 


-0.201 


-0.087 


0.083 


0.009 


Caha 


506.4 


23.8 


-0.207 


-0.085 


0.093 


0.011 


Caha 


510.5 


27.9 




-0.292 






WD 


511.4 


28.8 




-0.286 


0.062 




WD 


512.6 


29.9 


-0.150 


-0.084 


0.199 


-0.218 


SWO 


516.5 


33.9 






0.078 




WD 


518.4 


35.8 




-0.037 


0.082 




WD 


519.4 


36.8 


-0.125 


-0.038 


0.120 


-0.304 


LT 


520.4 


37.8 


-0.196 


-0.076 


0.076 




WD 


521.4 


38.8 




-0.078 


0.127 


-0.304 


LT 


524.4 


41.8 


-0.191 


-0.054 


0.074 




WD 



" JD - 2 453 000.0 Epoch in days with respect to the esti- 
mated B-band maximum JD 2 453 482.6 ± 0.5. " See Tabled 
for details. 




Figure S. B — R difference image of NGC 4070, constructed as 
described in the text. Dark shades in the figure correspond to 
red areas. The field of view is 2 X 2 arcmin^; north is up, east 
to the left. Crosshairs mark the explosion site of SN 2005bl and 
the centre of NGC 4070. 'B' denotes a bluish patch, probably a 
companion galaxy of NGC 4070 (see discussion in the text). 



iPatil et al.ll2007l for a more sophisticated picture). This is 
consistent with the negligible host-galaxy extinction found 
in SNe 1991bg, 1997cn and 1998de. However, recent merg- 
ers or interaction may alter this picture, as exemplified 
by the radio galaxy Centaurus A, the host of SN 1986G 
(|Phillips et al.l Il987l ). In this case, the SN lay behind a 
prominent dust lane, and consequently was strongly extin- 
guished. In NGC 4070 no such dust lane is visible even 
in deep images, but SN 2005bl showed signs of extinction 
within its host galaxy, the most prominent being a narrow 
interstellar NalD line in the spectra at the redshift of the 
host, with an equivalent width (EW) of 2.6 ± 0.3 A. 

Since elliptical galaxies lack Hll regions with strong 
emission lines, their surfac e colours are goo d tracers of the 
internal dust distribution l|Patil et al.ll2007l ). Hence, to in- 
vestigate the dust content in NGC 4070, we constructed a 
B — R image of the galaxy from the templates obtained with 
CAFOS on UT 2006 March 28, following largely the pre- 
scription of lPatil et al.l (|2007h : after the usual pre-reduction, 
the B- and 7?-band images were spatially aligned, and 
the sky background was subtracted. The images were then 
scaled to contain the same fiux inside an aperture of 45 
arcsec around the centre of NGC 4070, and subtracted one 
from the other. The difference image is shown in Fig.O Both 
the centre of the galaxy and the position of SN 2005bl are 
marked. Dark shades correspond to redder areas, indicative 
of either dust or an intrinsically redder stellar population. A 
red region is present to the immediate east and north-east 
of the nucleus, and another more extended but less opaque 
arc to its south-west. This asymmetric surface-colour distri- 
bution suggests that there probably is dust in NGC 4070. 
However, at the exact position of SN 2005bl no major blue or 
red structures can be discerned, so that the dust obscuring 
the SN is probably too locally confined to be resolved. 

Deep images of NGC 4070 reveal some deviation from 
a perfectly spherical or ellipsoidal shape (Fig. [Jjl. This is 
an indication of fairly recent interaction, either with the 
galaxy 2MASX J12040831+2023280 about 1.3 arcmin to its 
south-west (labelled 'A' in Fig. |4)|, or with a small knot 
about 0.8 arcmin to its south (labelled 'B', and most easily 
seen as irregularly-shaped bright patch in the lower part 
of Fig. |3]). Furthermore, in Fig. |4] a faint, broad bridge 
of luminous matter between NGC 4070 and its equally 
massive elliptical neighbour galaxy NGC 4066 can be de- 
tectedlfl The two galaxies are offset by about 3.74 arcmin, 
which corresponds to a projected distance of 114 kpc (for 
Hq = 72kms-iMpc"^). 

The exact amount of dust extinction towards SN 2005bl 
is quite difficult to determine, and it constitutes the main 
uncertainty in the calibration of the SN absolute magni- 
tudes. The Galactic comp onent is small, with a colour ex- 
cess E{B-V) Ri 0.03 mag l|Schlegel et al.lil998l '). However, as 
mentioned before, the contribution of dust in the host galaxy 
is significant. With EW(NalD)= 2.6 ± 0.3 A and applying 

E{B-V) = 0.16 X EW(NaiD), 

jTuratto. Benetti fc CappeUarol bOOSl ). we obtain E{B - 



^ The situation is so mewhat reminiscent of SN la 2005cf 
llPastorello et al.ll2007ah . with the difference that this SN was 
directly located in the tidal bridge. 
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Figure 4. Deep 7.7 X 9.0 arcmin-^ exposure of NGC 4070 and 
its neigiibour NGC 4066 obtained on UT 2006 March 28 with 
CAFOS. North is up, east to the left; brightness increases from 
blue over green to red. The distorted shape of NGC 4070 and 
the bridge of luminous matter connecting the two galaxies are 
discernible. 'A' and 'B' mark companion galaxies of NGC 4070. 



gous to that presented bv lElias-Rosa et al] (|2006l ) and based 
on extracting an extinction curve from a comparison with co- 
eval spectra of SNe 1991bg and 1999by, yields a host-galaxy 
colour excess of 0.22 ± 0.02 mag. Furthermore, no deviation 
of the total-to-selective-extinction parameter Rv from the 
canonical value of 3.1 as inferred from "standard" dust in 
the Milky Way can be discerned. 

For the rest of the discussion we assume E{B—V)host = 
0.17±0.08 mag. This estimate is based solely on the study of 
the SN colours (giving the strongest weight to the compar- 
ison with SN 1991bg), and ignores the larger colour-excess 
estimate from the interstellar NalD line, since the latter 
would result in too blue a colour and too high an absolute 
luminosity given SN 2005brs spectrophotometric similarity 
to SNe 1991bg and 1999by. This choice of E{B-V)host, to- 
gether with the foreground reddening of 0.03 mag, yields 
a total colour exces s E(B — V)totai = 0.20 ± 0.0 8 mag and, 
adopting a standard lCardelli. Clayton fc MathisI (Il989l ) red- 
dening law with Rv = 3.1, a total _B-band extinction along 
the line of sight of Ab = 0.82 ± 0.33 mag. 

With a recession velocity corrected to the CMB refer- 
ence frame of 7534 ± 33 kms"^ (Tabled, NGC 4070 is well 
within the Hubble flow. Adopting Hp = 72kms~^Mpc~^ 
iFreedman et al.|[200ll : ISpergel et al.ll2003l '). this corresponds 
to a distance of 104.6 Mpc and a kinematical distance mod- 
ulus fi = 35.10 mag, similar to that of the Coma Cluster. 
However, the latter is about 14°45' away (which, at the given 
distance, corresponds to 27 Mpc), excluding any physical as- 
sociation. Accounting for an uncertainty of 300 kms~^ aris- 
ing from a possible peculiar motion of NGC 4070, we obtain 
= 35.10 ± 0.09 mag. 



V^)host = 0.42 ± 0.05 mag. However, the lack of knowledge 
of the extinguishing material's gas-to-dust ratio makes this 
method very uncertain. Therefore, in the case of Type la 
SNe, alternative ways to determine the colour excess using 
the light and colour cur ves a r e usua lly preferred. Among the 
most well-known is thelLiral l| 19951 ) relation, which assumes 
a uniform B — V colour evolu tion of SNe la bet ween 30 and 
90 d after maximum (but see lWang et al.l [20071 for a recent 
warning about the use of this method). Apparently, this re- 
lation seems to hold also for underluminous SNe la, which 
are characterised by much redder B — V colours at maximum 
and soon thereafter. However, the decreasing quality of the 
SN 2005bl photometry after -1-30 d and the lack of B-band 
observations later than -1-42 d, result in a relatively large 
range of colour excesses consistent with the Lira relation, 
from E{B-V)i,ast « 0.12 to 0.25 mag. 

Since SN 2005bl is an underluminous SN la with 
Ami5(B)true = 1.93 (cf. Section [4T|, SNe 1991bg and 
1999by (Ami5(B)truo = 1.94 and 1.90, respectively; cf. Sec- 
tion [7]) are natural comparison objects. Despite their nu- 
merous similarities, the reddening-corrected colour curves 
of the two latter objects differ especially at early phases, 
SN 1991bg being ~0.15 mag redder in B - V than 
SN 1999by. This intrinsic colour difference directly prop- 
agates to an uncertainty in the inferred colour excess of 
SN 2005bl determined on the basis of such comparison. 
Matching the colours of SN 2005bl to SN 1991bg, we ob- 
tain E{B—V)hoBt ~ 0.13 ± 0.05 mag, whereas a comparison 
with SN 1999by yields £(B-"l/)host « 0.29 ± 0.05 mag. 

A study of the extinction law towards SN 2005bl, analo- 



4 PHOTOMETRIC EVOLUTION 



4.1 Filtered light curves 



The photometric observations of SN 2005bl are among the 
earliest ever obtained for a 91bg-like SN. In Fig. [5] we 
present the Bessell UBVRI and Sloan gz light curves, i.e. 
the data of Tables [3] and O and Contreras et al. (in prep.). 
5* -I- /f-corrections have been applied to the BVRI bands. 
Also shown are _BV7- templates constructed from SN 1991bg 
(Ha muv et all 1 199601. reddened by E{B - V) = 0.20 mag. 
These provide an excellent match to SN 2005bl. 

Compared to normal-luminosity SNe la the light curves 
of SN 2005bl are characterised by a fast rise to and de- 
cline from the light-curve peak, especially in the blue 
bands. Moreover, in the B band the settling to the ex- 
ponential tail (corresponding to the bend parameter of 
iPskovsklHl 19841 and the intersection parameter if as defined 
bv lHamuv et al.lll996d) occurs at a remarkably early epoch, 
only about 15 d after maximum, as compared to 25-38 d 
for intermediate and slow decliners (|Hamuv et al.lll996d ). 
Consequently, the decline from the peak to the onset of the 
radioactive tail is only ~ 1.9 mag, less than in most SNe la 
with shallower initial decline^ 

The exact value of the iPhiUipj ||1993D decline rate pa- 
rameter Ami5(_B) is difficult to measure directly because of 
a gap in the light curves around -1-15 d. A polynomial fit in- 
cluding the B-band data up to one month past maximum 
yields Ami5(_B) = 1.91, and we adopt a conservative error 
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Figure 5. UBVRI Bessell and gz Sloan light curves of SN 2005bl. 
The BVRI data (Table [Sjl are S- and ii'-corrected except for the 
latest phases, while the U, g and z data ( Ta bles \3\ \5 \ and Contreras 
et al. in prep.) are not. Hamuy et al.'s lll996cl) templates for 
SN 1991bg are shown for comparison (dotted lines). 



of 0.10 to account for the fact that the fit around +15 d is 
not very well constrained. After correction for the mitigat- 
ing effects of extinction (Phillips et al. 199!|) this turns into 
an actual decline rate Ami5(B)truo ~ 1.93 ±0.10. 

Besides the B band, differences with respect to normal- 
luminosity SNe la are most pronounced in the near-IR (/ 
and z bands), where - like in other 91bg-hke SNe - the sec- 
ondary light-curve maximum is absent, and the main maxi- 
mum is delayed with respect to that in B (JD 2 453 482.6 ± 
0.5) by a couple of days rather than advanced. In general, 
the instance of peak brightness seems to be the more delayed 
the redder the band is. T/-band maximum occurs 2.3 d after 
that in B, and those in R and / are delayed by 3.3 and 4.4 d, 
respectively (see Table [1]), whereas in U the light-curve peak 
precedes that in B by about 1.0 d. For the z band the sparse 
photometric coverage does not allow for an exact determina- 
tion of the time of maximum light, but it can be estimated 
to be similarly delayed as in R or I. 

The absolute peak magnitudes of SN 2005bl (Table [TJ, 
calculated adopting the distance and extinction estimates 
presented in Section O reveal that the SN is underlumi- 
nous by 1 to 2 mag in all filters compared to a canonical, 
Ami5(i?)truc ~ 1.1 SN la, the difference being most pro- 
nounced in the blue bands. 



4.2 Bolometric light curve 

The differences in luminosity can clearly be seen in Fig. |S1 
where the quasi-bolometric light curve of SN 2005bl is 
compared to those of the other underl uminous SNe la 
1991bg 1997cn, 1998de and 1999bv jFilippenko et"all 
1992bl: iLeibundgut et all Il993l: [ Turatto et al. 1 1 19961 . 1 19981 : 
Modiaz et al. 2001 : Garnavich et al.l 20041 ). and the in- 
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Figure 6. Quasi-bolometric light curves of SNe 2005bl, 1991bg, 
1997cn, 1998de, 1999by, 2004eo and 2005cf, obtained by integrat- 
ing the l/-through-/-band fluxes (for the adopted distance a nd 
extinction parameters see Table l8l and IPastorello et al.|[2007al lb[) . 
Error bars are shown for SN 2005bl only, and account for uncer- 
tainties in the photometric calibration, distance and extinction 
estimates. The Ami5(iJ)truo of the SNe is given in parentheses. 



termediate decliner SN la 2004eo (Ami5(_B)truc ~ 1.46, 
IPastorello et al ] |2007bh . which forms a bridge to canonical, 
norm al-luminosity SNe la su ch as 2005cf (Ami5(_B)truc = 
1.12, IPastorello et al. 2007a). The integrated optical light 



curves (see, e.g.. iNomoto. Filippenko fc Shigevamal I199GI ) 
were constructed in the following way: in a first step the 
?7-through-/ magnitudes were converted to monochromatic 
fluxes and the spectral energy distribution (SED) was in- 
terpolated linearly. The SED was then integrated over fre- 
quency, assuming zero flux at the integration limits, which 
are given by the blue edge of the U band and the red edge of 
the / band. Whenever no [/-band observations were avail- 
able, or the coverage of this band was incomplete, a correc- 
tion derived from SN 1999by was applied to the analogously 
constructed B-through-/ light curve. This method appeared 
more reliable than applying U -band corrections based on or- 
dinary SNe la. 

As Fig. [6] shows, 91bg-like SNe form a fairly homoge- 
neous group in terms of bolometric light-curve shape and 
luminosity. Both their light-curve width and their luminosi- 
ties distinguish them even from intermediate decliners such 
as SN 2004eo. SNe la with Ami5(B)truo ~ 1.50-1.85 (which 
are rare, cf. Section 17. ip would probably fall in the gap 
between SN 2004eo and the 91bg-like SNe. In Fig. [6] they 
are not included because their light-curve coverage is in- 
sufficient in some of the relevant bands. Within the 91bg- 
like group, SN 2005bl appears to be the brightest object 
(logLinax = 42.31 ± 0.11), but the differences are mostly 
within the error bars, which for SN 2005bl are dominated 
by the uncertainty in the host-galaxy extinction. 

In the only rapidly-declining S N la with ex t ended 
near-IR photometry SN 1999by (|Hoflich et al.l |2002| : 
iGarnavich et al.ll2004l fl. the JHK bands contribute ~ 25 % 
to the total bolometric flux around maximum, and ~ 36 % 

9 SN 1991bg has been observed in JHK bv lPorter et al.l lll992l ') 
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Figure 7. Time-evolution oi the U - B , B ~V , V - R and V -I 
colour indices of SN 2005bl. The c u rves have been reddening- 
corrected adopting a ICardelli et alj lll989l ) extinction law with 
E{B -V) = 0.20 mag and Ry = 3.1. 

after one month. The corresponding numbers measured for 
SN 2004eo are ~ 16 % and ~ 37 %, respectively. This means 
that at early phases the near IR gives a larger contribution 
in the fast decliners than it does in intermediate and slow 
decliners, but that by -fSOd the difference has vanished. 
Wavelength regions other than the optical and near-IR ap- 
pear to con tribute very li t tle to the total bolometric flux 
in SN e la (jSuntzcfl 1996; IContardo. Leibundgut k, Vaccal 
Therefore, neglecting these regimes should not cause 
a significant underestimate of the true bolometric luminos- 
ity. 

4.3 Colour evolution 

Fig. [7] presents the time-evolution of the U — B, B — V , 
V - R &nAV - I colours of SN 2005bl. Th e good agree- 
ment with the colour curves of SN 1991bg (|Hamuv et al.1 
Il996d ) also shown in the figure is evident. Except for the 
late V — R points, all colour indices are > throughout the 
investigated period, indicating that SN 2005bl - like other 
rapidly-declining SNe la - was a rather red event. In particu- 
lar it was much redder at early phases than ordinary SNe la, 
which are characterised hy a. B — V between 0.0 and —0.1 at 
maximum light. Remarkably, the evolution of SN 2005brs 
various colour indices with time is very similar. From our 
first observations at —6 d on, all colours become monotoni- 
cally redder until at least ~ 10 d after maximum, and almost 
simuhaneously exhibit a red peak between -1-12 and -1-17 d, 
followed by a monotonic bluening which probably lasts be- 
yond the end of our photometric coverage around -1-40 d. 

(five epochs), but only the three ep ochs close to maximum light 
have been calibrated and published llKrisciunas et ahllioQ^) 
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Figure 8. Time sequence of SN 2005bl spectra in the observer 
frame. The phases reported next to each spectrum are with re- 
spect to _B-band maximum. The —3d (blue part) and -|-33d spec- 
tra have been smoothed using kernel sizes of 600kms~^. Gaps in 
the TNG spectra indicate the positions of the strongest night-sky 
residuals, which have been cut in order to avoid confusion with 
true SN or host-galaxy features (see Section 12.21 and Table |6ll. 

Such a high degree of similarity between different colour in- 
dices is not encountered in normal-luminosity SNe la (cf. 
Section rn^. 



5 SPECTROSCOPIC EVOLUTION 
5.1 Spectra of SN 2005bl 

The spectra of SN 2005bl presented in Fig. [S] cover the in- 
terval from 6d before to 33 d after B-band maximum hght, 
i.e., the photospheric phase and the early transition stages 
towards the nebular phase. During this entire period, the 
SN follows the evolution of other SNe la, in particular of 
the 91bg-like subclass. The pre-maximum spectra have a 
blue continuum, with characteristic P-Cygni lines of Sill, 
Sll, Call and Mgll superimposed. With respect to normal- 
luminosity SNe la, additional strong Till absorptions are 
visible, and the Si ll A5972 and O I A7774 absorptions are 
more pronounced. Shortly after maximum hght the blue fiux 
decreases significantly, and also the S ll and Si ll A5972 lines 
fade rapidly, and are no longer detectable two weeks after 
maximum. Si ll A6355 is, like most other photospheric lines, 
somewhat more persistent, and can still be discerned beyond 
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Figure 9. Comparison of prc-maximum spectra of underlumi- 
nous SNe la. The spectra were reddening-corrected according 
to the B — V colour excesses reported in Section [7] The — 6d 
spectrum and the blue part of the combined —3d spectrum of 
SN 2005bl were boxcar smoothed using kernel sizes of 2300 and 
3400 km s~^, respectively. 

+20(1. Nal D shows a trend opposite to that of Si and S, 
first being just visible as a shoulder in the blue wing of the 
Sill A5972 line on day +5, but evolving to a distinct absorp- 
tion feature by day -1-20. In the same period, Fell emission 
lines start to dominate the spectrum. 

5.2 Spectroscopic comparison with other 
underluminous SNe la 

In Figs. |9] to [11] the spectra of SN 2005bl are compared 
with those of other underluminous SNe la at epochs of 
about — 5d, -1-4 d and -|-30d, respectively. For illustration, 
a comparison with the spectra of t he intermediate decliner 
SN la 2004eo (Ami5(-B)truo = 1.46. iPastorello et ahllioOTbl ) 
is made in Fig. 1121 The figures show that 91bg-like SNe 
form a relatively homogeneous spectroscopic subclass, dis- 
tinct from SNe la with normal luminosity and even from 
those with Ami5(_B)truo close to 1.5. 

Already at —5 d (Fig. [9} the underluminous SNe la 
show evident Ti features between 4000 and 4400 A, which 
are absent in SN 2004eo (Fig. I12p where this region is dom- 
inated by Si and Mg lines. At the same time Fig. [12] sug- 
gests that at — 6d the continuum temperatures are not too 
different, indicating that the strength of the Ti features in 
91bg-like SNe may not be a pure temperature effect but 
require a trul y larger Ti abundan ce (see also Section 16. 2|) . 
In SN 1986G (Phillips et al."l987") the Ti lines are less pro- 
nounced than in S Ne 1998dc (Mathcson et al. 20071), 1999by 
iGarnavich et al] l2004i ) and 2005bl, emphasising its tran- 
sitional character between "normal" and strictly 91bg-like 
SNe la. Also, the absolute depth of the Sin A6355 line in 
SN 1986G resembles more that of SN 2004eo than those 
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Figure 10. The same as Fig. [9] but some days after maximum 
light in B. 
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Figure 11. The same as Fig. [9| but 4-5 weeks after maximum 
light in B. The spectra of SNe 1997cn and 2005bl were smoothed 
using kernel sizes of 600kms~^. 

of SNe 1999by and 2005bl. The Sill lines and most other 
features of SN 1998de at day —5 are broader and at bluer 
wavelength than in the other SNe of the sample. The spectra 
of SNe 2005bl and 1999by are essentially identical, the only 
exception being the flux depletion in SN 2005bl redwards 
of the Sin A6355 line at —6 and — 5d, which we attribute 
to Cll A6580 (see also Section \T^ . Ol A7774 is particu- 
larly pronounced in all underluminous SNe la for which the 
spectral region is covered. 
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Figure 12. Spectroscopic com parison of SN 2005bl an d tlie in- 
termediate decliner SN 2004eo ('Pastore llo et al.ll2007bl) at three 
different epochs. The spectra of SN 2004eo have been dereddened 
for E{B -V) = 0.109 mag; the -6d spectrum of SN 2005bl was 
smoothed using a kernel of 2300 km s~^. 
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Figure 13. Expansion velocities of Sll A5640, Sill A6355, Ol 
A7774 and the Can NIR- triplet as measured from the minima of 
the P-Cygni line profiles in SN 2005bl. 



at late epochs llFilippenko et al" l ll992bl : lTuratto et al. I ll996l : 
iMazzali et al■lll997^ ■ is also present in SN 2005bl. 



By a few days after maximum light (Fig. [10} the spectra 
have evolved significantly with the continuum being much 
redder now, but the degree of homogeneity is still remark- 
ably high. The Ti troughs in the blue are now fully de- 
veloped, showing the characteristic flat bottom that dis- 
tinguishes 91bg-like from other SNe la (Fig. [H]). The W- 
shaped S ll lines around 5500 A are comparatively weak. 
Again SN 1986G takes an intermediate position, with the 
characteristic properties of the underluminous class being 
less pronounced than in the other objects shown in Fig. 1101 

The spectra taken 4-5 weeks past maximum (Fig. Illf) 
reveal that at those phases the transition to the nebu- 
lar phase has already started. Emission lines of Fe-group 
elements are visible, but the pseudo-continuum has not 
yet vanished. Call, Ol and the Ti trough are still promi- 
nent in absorption, while S ll lines cannot be identified any 
longer. A characteristic feature of many 91bg-like SNe la 
at those epochs is the remarkably narrow NalD absortion 
near 5700 A. Especially in SNe 1991bg (Filippcn ko et al.1 
ll992bl:lLeibundgut et ahll 19931 : iTuratto et al. 19961 1997cn 
( Turatto et al.l Il998l ') and 1999by this line is very deep 
and distinct, while it is less pronou nced in SNe 2005b l 
and, in particular, 1986G and 1998de ( Modiaz et al.ll200lh . 
where also th e narr ow core is absent. As already noted by 
iModiaz et al] (1200 ll ). SN 1998de deviates significantly from 
other 91bg-like SNe between 6800 and 7700 A, exhibiting a 
relatively smooth continuum without strong O I A7774 fea- 
ture. 

Unfortunately, no late-time spectra of SN 2005bl were 
obtained, so that it cannot be verified whether the distinct, 
broad nebular emission feature near 7300 A, which charac- 
terises SN 1991bg and distinguishes it from ordinary SNe la 



5.3 Ejecta velocities 

Although ejecta velocities inferred from the blueshift of the 
absorption minima of P-Cygni lines suffer from a number 
of uncertainties (such as the typically strong line blending 
in SN la spectra), they do provide important information 
about the kinetic energy of the ejecta, and can be used as an 
observationally accessible parameter for comparison studies 
between different SNe. Furthermore, the range of velocities 
encompassed by different elements provides insight into the 
chemical stratification of the ejecta, and hence into nucle- 
osynthesis conditions. This fact has recently be en made use 
of in the Zorro diagnostics (|Mazzali et al. I l2007ri . In SNe la, 
especially the Si ll A6355 has proved to be a suitable veloc- 
ity indicator, since the line is usually well pronounced, fairly 
unblended and well visible over a relatively long period. 

The Si line velocity in SN 2005bl (Fig. [T3j evolves from 
- 11000 kms~i at -6d to ~ 6600 kms"^ at -|-33d. At 
maximum light, it is about 10 000 kms~^, which is low for a 
SN la. The velocity evolution of S ll A5640 resembles that of 
Sill, but the temporal decrease is steeper, and the absolute 
values are systematically lower by 500-2500 kms~^, mean- 
ing that the line predominantly forms in deeper layers of the 
ejecta (note that the line is visible only until -1-5 d). The Call 
IR triplet, on the other hand, has higher velocities than Si ll, 
and the difference increases from about 1000 kms~^ at — 6d 
to 3500 kms~^ at +33 d. This relatively shallow decrease in 
velocity means that Can lines mainly form well above the 
photosphere. Finally, O I A7774 also exhibits a fairly con- 
stant velocity of 9000-10 000 kms~^ during the entire pe- 
riod of our observations, but this result is not very robust as 
the feature may be substantially blended with Mgll A7890. 
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6 SPECTRAL MODELLING 



A ID Monte-Carlo spectrum synthesis code (|Mazzalill2000l ) 
was used to simulate the radiation-transport processes in the 
expanding ejecta of SN 2005bl. Basic assumptions include 
spherical symmetry, a Chandrasekhar-mass explosion, and 
an under lying density profile adopted from the W7 explosion 
model of lNomoto. Thielemann fc Yokoil (jl984l ). No attempt 
was made to obtain better results altering these ingredients. 
Synthetic fits to four early-time spectra of SN 2005bl were 
obtained and are presented here, after a short introduction 
to the concept of the code and the underlying model. 



6.1 Concept of the radiative transfer code 

Only a brief outline is given here; f o r detailed de scriptions of 
the code seelAbbott fc Lucy! (ll985l '). lMazzaH fc Lucv (199^. 
ILucvI (|l999l ) and iMazzalil \2Qm . where the basic develop- 
ments are documented. 

The radiative transfer is performed above the photo- 
sphere, which is located at an adjustable radius rph. Energy 
deposition from radioactivity is assumed to occur below the 
photosphere, from which a blackbody continuum is thought 
to be emitted. In the "atmosphere" above, radiative equi- 
librium is supposed to hold. Photons interact with lines and 
scatter on free electrons, but no continuum formation is as- 
sumed (Schuster-Schwarzschild approximation). The under- 
lying W7 density profile is scaled to match the epoch of each 
spectrum, assuming homologous expansion, i.e. r = vt for 
each particle^ Element abundances are assumed to be ho- 
mogeneous inside the envelope. They can be freely adjusted 
in order to match a given observed spectrum. 

In our Monte Carlo calculation, radiation packets are 
followed from their emission at the photosphere through 
their interaction history, until they either escape or are re- 
absorbed at the photosphere. The code takes into account 
scattering processes (on atoms / ions and electrons) as well 
as photon branching (absorption in an atomic line and sub- 
sequent reemission in another). Line optical depths are cal- 
culated in the Sobolev approximation, applicable for fast- 
expanding atmospheres. 

lonisation and excitation conditions are calculated from 
the radiation-field statistics. This is done using approximate 
NLTE formulae (see references above), which in principle 
employ LTE radiative rates, but additionally take into ac- 
count the dilution of the radiation field. Collision processes 
are neglected. A radiation temperature Tr, determining ra- 
diative rates, is calcula ted in each zone of the envelope. As in 
iMazzaU Lucvl l| 19931 ) it is chosen such that the mean fre- 
quency of the radiation field inside the zone matches that of 
a blackbody at Tr . The matter state and the radiation field 
are iterated until sufficient convergence is achieved. Within 
this process, the temperature of the inner-boundary black- 
body spectrum is adjusted in order to match the bolometric 
SN luminosity. 

Finally, the emitted spectrum is recalculated solving the 
formal integral, employing source functions obtained from 
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This is almost exact for our purposes, as the expan sion is 
homologous already ~ 10 s after explosion onset (see e.g. iRopkel 
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Figure 14. Synthetic fits to the SN 2005bl spectra at days —6.0, 
—5.0, —3.0 and +4.8. The observed spectra were smoothed with 
kernel sizes of 2300 kms"! (-6.0 d) and 3400 kms"! (-5.0 d 
and — 3.0d-blue part) for presentation. The values of relevant fit- 
parameters are summarised in Table [7] Atop the — 6d spectrum 
an identification of the most important lines is given. 



the packet statistics. This yields smooth spectra at relatively 
low packet numbers. 



6.2 Physical parameters and the chemical 
composition of the models 

The basic parameters and composition inferred from the 
best-fitting synthetic spectra are reported in Table [7] In 
Fig. 1 141 these models are compared to the observed spectra, 
and line identifications for the —6 d spectrum are shown. 
A rise time to the B-band maximum of 17.0 d is assumed, 
slightly less than the fiducial value of ~ 19.5 d for norma l- 
luminosity SNe la (Ricss ct al. 1999; IConlev et al] 120061 ). 
This choice is discussed in Section [7.21 

While the overall shape of the continua is nicely 
matched and most features are well reproduced, the main 
shortcomings of the models are a flux excess in the emis- 
sion component of Sill A6355 at early times, and, more 
worrysome, a mismatch in the position of some lines. 
Additionally, many synthetic features are too broad and 
strongly blended, as they were in models for SN 1991bg 
l|Mazzali et al.l Il997l ). No substantial improvement is ob- 
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Table 7. Physical parameters of SN 2005bl inferred from synthetic spectr a, assuming a rise time of 17 d in B. The mass fractions of 
selected elements are also recorded, and solar photospheric abundances ( As plund. Grevesse &: Sauvalll2005l ) are shown for comparison. 
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bolometric luminosity [ergs"^] 
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photospheric velocity [kms^ ^] 
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X(Mg) 
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X(Si) 


0.025 




0.037 
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tained reducing the photospheric velocity Wph, which indi- 
cates that the W7 density profile may not be perfectly suit- 
able for 91bg-like SNe, and that a steeper density gradient 
may be required. With Uph = 7500 kms"^ 6 d prior to maxi- 
mum light, the inferred photospheric velocity is significantly 
lower than in normal SNe la (~ 10 000 to 11000 kms"^ at 
comparable epochs). By day +5, Wph further decreases by 
about 1500 kms~^. Singly and doubly ionised species dom- 
inate the ejecta. Especially the singly ionised species leave 
strong imprints on the spectra, as can be seen from the line 
identification in Fig. 1141 In all spectra modelled here the 
heavy-element content is comparatively low, but Table [7] 
shows that it increases with time, indicating composition 
stratification in the ejectaQ At the higher velocities the 
ejecta are almost entirely made up of unburned material, 
whereas NSE elements are essentially absent. In particular, 
no Ni and Co, and only very little Fe are included in the —6 
and — 5d models. After maximum light, IMEs and Fe-group 
elements become more abundant, but still less so th an in 
normal- luminosity SNe la (cf., e.g.. IStehle et aL I I2OO5I ). In- 
terestingly, Ti and Cr, whose abundances can be well con- 
strained from the depth of the characteristic trough between 
4000 and 4400 A, appear to be more abundant than Ni, Co 
and Fe, yet another indication for explosion conditions which 
disfavour burning to NSE. 



7 DISCUSSION 



In a number of 


works (e.g. 'FilioDenko et al.| Il992b; 


Leibundgut et al.l 1 19931: iRuiz-LaDuontc et al.l 1993; 


Turatto et al. 1996 


: iMazzali et al.l Il997l: iTuratto et al. 


19981: iModiaz et all 


2OOII: iGarnavich et al.(l2004l) 91bg-like 



SNe la have been studied and found to be different from 
normally luminous SNe la in several respects. Besides their 
lower luminosity, they have rapidly declining light curves 
which do not show the characteristic secondary maximum in 



the near-IR, cooler spectra with significant Ti 11 absorption 
troughs, lo w ejecta velocities an d large values of the 7?,(Si) 
parameter (iNu gent et al] Il995l ). Nevertheless, the light- 
curve width and peak lumi nosity of underluminou s SNe la 
seem to obey a correlation (|Garnavich et al]|2004h . but not 
the same one as more ordinary SNe la. Also, 91bg-like SNe 
seem to fit smoothly into the Zorro plot l|Mazzali et al.l 
12007.) . where they represent the extremely ^^Ni-poor, 
IME-rich end of the SN la distribution (but see Section [7.41 
for a revised picture). With the increasing number of 
well-observed objects of this class, it is interesting to revisit 
some of the aspects mentioned above and to investigate in 
more detail the degree of homogeneity, or rather diversity, 
that these SNe exhibit. In Table [S] elementary information 
on our comparison sample of rapidly-declining SNe la 
is collected; selection criteria were a reasonably dense 
light-curve coverage starting at least around maximum 
light, and Ami5(i3)truc > 1.5. 



7.1 Photometric behaviour of underluminous 
SNe la 

Here we focus on the photometric properties of underlumi- 
nous SNe la, compare their light-curve shapes and colour 
indices, determine their peak absolute magnitudes, and in- 
vestigate their decline-rate vs. luminosity relationship. 



Light- and colour-curve morphology 

Fig. [15] compares the B, V and I light curves of all the 
objects in Table |8] rescaled to coincide at maximum light. 

The SNe with Amis (i?)truc ~ 1-9 (red symbols in 
Fig. I15|l , including SN 2005bl, show a single-peaked broad 
7-band light curve, whose maximum is delayed by a few 
days with respect to that in B. In contrast, the SNe 
with Amis (B)true ~ 1.5 to 1.75 (green symbols) clearly 



For this reason it w ould be necessary to work with the strati- 
fied version of the code llStehle et al .120051) if accurate abundances 
were required. This reaches beyond the scope of this paper. There- 
fore, the composition reported in Table[7]refiects the overall trend, 
but the absolute numbers should be taken with caution. 



For SN 1986G the Cepheid measurement does not refer to 
the actual host galaxy Cent aurus A, but to NGC 5253, which is a 
member of the same group CF errarese et al.l200Clt) . It was adopted 
here for the lack of more reliable distance estimates. 
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Table 8. Comparison sample of SNe la with Ami5(B)truc > 1-5, ordered by increasing Ami5(_B)true- 



SN 


Ami5(B)truo 


Ms, max 




7e(Si)max 




E{B-V)^ 








morph.'* 


references 


1990af 


1.57(0.05) 


-18.96(0 


25) 


-18.98(0 


19) 








0.07(0.06) 


36 


59(0 


05) 


SBO 


H96a, H96b, P99 


ZUUIJUK 


1 c^v/fi ^n^ 

i.o ( (u.uy ) 


-18.84(0 


18) 


-18.83(0 


17) 








n (\'7(r\ nQ\ 
Kj.yj ( yJ.yjo ) 


o4 


1 c /n 


13) 




JUo, VoiNrji 




1 f^Q/Tl 0^^^ 
± . uy ^^u . uo j 


-18.60(0 


30) 


-18.73(0 


28) 








yj.uuyyj.uo j 


32 


82(0 


25) 


E2 


Tflfi VSNF-T 

JUU, V OiN H/ J- 


1992bo 


1.69(0.05) 


-18.61(0 


17) 


-18.59(0 


15) 








0.03(0.03) 


34 


34(0 


12) 


SBOpcc 


H96a, H96b, P99 


1993H 


1.70(0.10) 


-18.57(0 


25) 


-18.73(0 


20) 


0.50(0 


05) 


73(8) 


0.12(0.06) 


35 


07(0 


09) 


SBab(rs) 


H96a, H96b, P99 


1986G 


1.81(0.07) 


-17.76(0 


32) 


-18.03(0 


26) 


0.60(0 


04) 


68(4) 


0.78(0.07) 
Rv=2.4<= 


27 


61(0 


11) 


SOpcc 


H87, P87, P99, 
FOO, B05 


1998bp 


1.83(0.06) 


-17.73(0 


25) 


-18.08(0 


24) 








0.08(0.03) 


33 


13(0 


22) 


E 


G04, J06 


1997cn 


1.88(0.10) 


-17.17(0 


20) 


-17.79(0 


19) 


0.68(0 


07) 


75(9) 


0.03(0.03) 


34 


25(0 


13) 


E 


T98, B05, J06 


1999by 


1.90(0.05) 


-17.17(0 


26) 


-17.65(0 


25) 


0.69(0 


05) 


97(4) 


0.02(0.03) 


30 


75(0 


23) 


Sb 


TOO, HOI, vol, 
G04, VSNET 


2005bl 


1.93(0.10) 


-17.24(0 


34) 


-17.85(0 


27) 


0.63(0 


06) 


99(9) 


0.20(0.08) 


35 


10(0 


09) 


E 


this work 


1991bg 


1.94(0.10) 


-16.85(0 


34) 


-17.56(0 


29) 


0.66(0 


05) 


106(5) 


0.08(0.06) 


31 


28(0 


20) 


El 


F92, L93, T96, 
P99, FOO, TOl 


1999da 


1.95(0.10) 


-16.98(0 


22) 


-17.65(0 


20) 








0.06(0.03) 


33 


58(0 


17) 


SAO 


KOI 


1998de 


1.95(0.09) 


-16.74(0 


19) 


-17.43(0 


17) 


0.69(0 


04) 


146(3) 


0.06(0.03) 


34 


06(0 


14) 


SO 


MOl, J06 



" Post-maximum decrease of the Sin A6355 velocity in kms ^ d ^, see lBenetti et aP l|2005h ans Section 17.31 Total (Galactic + host 
galaxy) colour excess. " Distance modulus from Cepheids (SNe 1986G and 1999byO SBF and PNLF (SN 1991bg) or the host galaxy 
recession velocity with respect to the CMB rest frame (NED) assuming Hq = 72kms~^Mpc~^ (other SNe). For the latter, an 
uncertainty of 300km s~^ has been adopt ed to account for the galaxies' peculiar motions. Host-galaxy morphology (from LEDA). 
" Spectropolarimetry l lHough et al.|[l987^ indicates Ry = 2.4 ± 0.13 for the dust in Centaurus A, the host of SN 1986G. We adopted 
this value and propagated the assigned error to the uncertainty in the absolute magnitudes. 



H9 6a =lHamuv et al.l l ll996al ') ; 
iHough et al.l lll987l'): P87 



H96b = iHamuy et al.l l ll996bl) ; P99 



H9bb =|JlamuY et al.l 
= |Phillips et al.l il987t); 



= iPhillips et al.l l ll999l) : J06 



FOO 



_ iJha et al.l ||2006|) : H87 : 
iBenetti et al. !|2005)); G04 



■ ■ ■ ■ , OOl); B05 

iGarnavich et al. " il200i) ; T98 =lTuratto et al.l lll998l) : TOO =lTbt h fc Szabo (2000'); HOI = iIowell ct al., (.20011'): VOl 
[•yinko et al.' (2001'); F92 =|Filippcnko ct al.' ('l992b'); L93 = iLeibundeut et al. (1993); T96 = lTuratto et all (Il996h ; TOl 
Irbnrv ct al. (200l] ); KOI = ,Krisciunas ct al.. C2001i) ; MOl = lModiaz et al. (2001); 'VSNET = "Variable Star Network 



have double-peaked 7-band light curves, and, as in normal- 
luminosity S Ne la, the first J-band peak precedes tha t of the 
B band (e.g. iHamuv et aLlll996d : lLeibundgutll2000l ) These 
observ ations are in agreement with the results of iKasenI 
(|2006l ). who computed synthetic light curves with the SE- 
DONA code, varying the mass of ^^Ni. His prediction is an 
ever smaller temporal offset of the first and secondary I- 
band maxima with decreasing SN luminosity, and the merg- 
ing to a single broad peak for the most underluminous ob- 
jects. The _B-band light curves of the two groups (red and 
green symbols in Fig. I15|) also exhibit noticeable differences 
in shape, probably influenced by the evolution of the Ti 
trough. 91bg-like SNe settle to the radioactive tail earlier, 
only ~ 15d after maximum light. Consequently, the drop in 
magnitude from the peak to the tail is smaller for these, 
although their initial decline is steeper. Thus, the B light 
curves of 91bg-like and "normal" SNe la cannot be trans- 
formed into each other by employing a simple strech factor. 

The only objects in the sample with 1.75 ^ 
Ami5(B)truo < 1.85 are SNe 1986G and 1998bp (yellow 
symbols). The light-curve morphology of these SNe seems to 
provide a link between the formerly defined groups. While 
the I band is still single-peaked with maybe a hint of a 
double peak, the luminosity drop from maximum to the ra- 
dioactive tail in the B band resembles that of SNe with 
Ami5(_B)truc < 1.75. The existence of objects with such in- 
termediate properties might support the idea of a common 
progenitor and explosion scenario for all SNe la. However, 
SN statistics tell us that objects with Ami5(i3)truc ~ 1-8 
are intrinsically rare. 



Fig. 1 16 1 shows the reddening-corrected B — V and V — I 
colour curves of the same SNe as above. The differences 
in B — V a fter +30 d are marginal, and all SNe obey the 
iLiral (|l995l l relation within the uncertainties (note, how- 
ever, that for the more strongly-reddened SNe 1986G and 
2005bl the Lira relation was directly or indirectly used to 
infer the true extinction along the line of sight, making this 
a circular argument). On the contrary, the B — V colour evo- 
lution around maximum light shows remarkable differences, 
the 91bg-like SNe la having a much redder colour at max- 
imum (0.4 to 0.7 mag, compared to ~0.0 mag for SNe la 
with Ami5(_B)truc < 1.75). Also, the peak B — V colour is 
reached earlier in underluminous SNe, and is redder by ~ 0.3 
mag. 



\n V ~ I the differences among the objects of Table [8] 
are even more pronounced. The V — I colour index of SNe la 
with 1.50 < Ami5(_B)truo < 1-75 decreases from about 10 d 
before to 5-10 d after B-band maximum just like in all in- 
termediate or slow decliners, then increases steeply, levelling 
off at ~ 20 d and again decreasing slowly thereafter. 91bg- 
like SNe, in contrast, do not show the initial bluening, but 
become redder from the very first available observations on- 
wards. This behaviour continues until ~ 15 d past maximum, 
when the curves flatten and subsequently a soft bluening sets 
in. Evidently, the origin of the different V — I evolution lies 
in the different delay of the main maximum and the presence 
or absence of the secondary maximum in the I band. 
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Figure 15. BVI light curves of the rapidly-declining SNe la of 
Table |8] normalised to their peak magnitudes. The Ami5(B)truo 
of each SN is given in parentheses. Different symbol colours rep- 
resent different decline rates (green: 1.50 < Ami5(i?)truc < 1.75; 
yellow: 1.75 < Ami5(i?)truc < 1.85); red: Ami5(i?)truc > 
1.85). Light-curve t emplates of SN 1992al {Ami5(B)truc = 1.11, 
iHamuv et"al]|l996d) are shown for comparison. 



Absolute magnitudes & Phillips relation 

Table |8] shows that the peak absolute magnitudes of the 
SNe with Anii5(B)true > 1-5 spa n a wide ran g e, from 
-16.74 to -1 8. 96 in the B band. AslHamtiv et all (|l996al ). 
iPhillips etahl ([l99i) and iGarnavich et alj (|20oi ) pointed 
out, also for these SNe the peak magnitudes correlate 
with Ami5(B)truo, but the dependence is steeper than for 
more slowly declining SNe la. Based only on SNe with 
Ami5(B)true ^ 1.69, our best linear fits (Fig. (TT]) are given 

by 



Ms.pcak = -18.54 5 log(Ho/72) + 6.83 (Amis (B)truc - 1-7) 

Mv.pcak = -18.61 + 51og(//o/72) + 4.33 (Ami5(B)truo - 1.7). 

The slopes of 6.83 ± 0.32 and 4.33 ± 0.31 for the B 
an d V bands are i n exce llent agreement with the results 
of 'Garnavich ct al.' ('2004'), and much steeper than what 
[Hamuy et al. (1996a ) find for SNe with Ami5(B)truc < 1.69 
(0.78 ±0.17 and 0.71 ±0.14 for B and V, respectively). Also 
the quadratic relation derived by IPhillips et al.l ()l999l l for 
SNe with Ami5(B)true < 1-70 provides a poor fit to the 
peak magnitudes of fast decliners. 




10 20 30 40 

epoch with respect to B maximum 

Figure 16. Evolution of the B — V (top panel) and V — I (bot- 
tom panel) colour indices of the rapidly-declining SN la sample 
of Table[8]plus SN 1992al. The symbols and colour coding are the 
same as in Fig. 1151 The dash-dotted line in the lower panel repre- 
sents the[Lira (1995) relation, i.e. the uniform B — V colour that 
all SNe la are supposed to exhibit between day -f 30 and -(-90. 
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Ann,5(B)„„, 

Figure 17. Peak absolute B- and V-band magnitudes of the SN 
sample of Table |8] The dotted lines are the best linear fits to 
the data with Ami5(B)truc ^ 1.69, characterised by slopes of 
6.83 ± 0.32 and 4.33 ± 0.31 for B and V, respectively The solid 
lines show the quadratic relation obtained bv iPhillips et al., (,199"9l l 
for SNe la with Ami5(i3)truo < 1.70. 
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7.2 The rise time of underluminous SNe la 

Rise times of SNe la are important to constrain possible ex- 
plosion models. In our synthetic spectra the rise time tr de- 
termines the actual (scaled) density profile at a given epoch 
(cf. Section [6. If) , and hence the matter density of the line- 
forming region. Here we report our attempts to estimate the 
rise time of 91bg-like SNe, in particular of SN 2005bl. 



Photometric rise-time determination 

Efforts have been made to measure SN-Ia ris e times di- 
rectly from a fit to their e a rly-time light cu rves l|Riess et al.l 
1 19991 : IConlev et all l2006l : IStrovinkl l2007h . Depending on 
measurement details, these studies yield rise times of 17 
to 20 d for fiducial Ami5(B)truc = 1.1 SNe la. Further- 
more, there seems to be a negative correlation between 
and Ami5(i?)true, faster decliners having shorter rise times. 
However, all these studies are based on samples of "normal" 
SNe la with Ami5(_B)truc ^ 1-5, and it may be doubted 
whether the inferred trends can be extrapolated to 91bg- 
like SNe with Ami5(B)truo « 1-9. 

What makes these direct measurements so demanding 
is the need for high-quali ty early-tim e photometry, prefer- 
entially before day —10 (|Riess et al.l 1999), which to date 
is not available for any 91bg-like SN, including 2005bl. In 
order to constrain the rise times of this class of objects, 
it may be useful to consider the two 91bg-like SNe la with 
the earliest photometric data, SNe 1998de and 1999by. Both 
have a filtered light-curve coverage starting about one week 
before B-band maximum, complemented by earlier unfil- 
tered measurements and some deep detection limits shortly 
before. For SN 1998de the earliest detection showing just 
"a hint of the supernova" was on unfiltered CCD frames 
taken 12.3 d before B-band maximum (preceded by a non- 
det ection down to a lim iting magnitude of 19.0 five days ear- 
lier, iModi az et al.ll 19981 ) , thus providing at least a lower limit 
for the rise time. For SN 1999by, the earliest detection on 
unfiltered CCD frames dates back to day —11. 4, showing the 
SN a t only 5 ± 1 % of the peak ii-Iuminosity (jArbour et al.l 
1 19991 ). Further detection limits constrain the SN luminosity 
to be less than ~ 1.6 % of the peak value at —12.4 d, and less 
than ~0.3% at — 15.1 d. Applying the method described by 
iRiess et al.l (| 19991 ) to these data, but using unfiltered and R- 
band points between —11.4 and — 4.0d for the fit (Fig. ll8p . a 
rise time to B-band maximum of 13.9ti i would be inferred 
for SN 1999by, the errors being the 3-a" confidence levels of 
the fit combined with the uncertainty in determi ning the 
epoch of B-band maximum (|Garnavich et al.ll2004l ). 

In this approach we disregarded a possible (and theo- 
retically not unexpected) deviation from the early L oc t^- 
behaviour which is the basis of the Riess et al. analysis. If 
the synthesised ^^Ni is confined to the inner ejecta, which - 
at least for 91bg-like objects - is supported by the lack of Ni 
in the -6 and -5d spectra of SN 2005bl (cf. Section [O)) . 
it should require some time for the photons to reach the 
photosphere and be released. This could result in an initial 
post-explosion phase with only very little brightening of the 
SN, before the observed steep rise of the light curve sets in. 
For this reason, the inferred value of 13.9 d might be con- 
sidered a lower limit to the actual rise time of SN 1999by. 
Given the similarity in terms of light-curve shape, this re- 
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Figure 18. R-band / unfiltered early-time light curve of 
SN 1999by, plotted as L^/^ vs. t - ts^^^, and linear fit (dot- 
ted line) to the data up to day —4 (between the two thin vertical 
lines). Assuming L oct^ , the instance of explosion would be given 
by the intersection of the fit with L^/-^ = at — 13.9d. The grey- 
shaded region marks the 3-(t confidence bands of the linear fit. 
Caveats are discussed in the text. 

suit should provide a good indication also for the rise time 
of SN 2005bl. 

The rise time addressed by spectral modelling 

The rise time constitutes an important input parameter 
for our spectrum synthesis calculations. Since the colour 
of an observed spectrum has to be reproduced in a corre- 
sponding model, the photospheric temperature, which can 
be crudely approximated through the Stefan-Boltzmann-law 
L oc Upjjt^T'', is confined to a limited range. Assuming a 
shorter rise time, the photospheric velocity thus usually has 
to be increased in order to keep the photospheric tempera- 
ture constant. The line velocities in the synthetic spectrum 
then also tend to be larger. 

Given the difficulties in observational rise-time deter- 
mination, we created eight models with different rise times 
(12.5-25.0 d) for the earliest available spectrum of SN 2005bl 
(— 6d), where a change in tr has the largest relative effect. 
Excluding models largely deviating from the observed spec- 
trum yielded a range of acceptable rise times. For the mod- 
els with tr < 14.0 d, we had to use photospheric velocities of 
«ph > 10 200 kms~^, which are incompatible with observed 
and measured (Fig. I13p O I and S ll expansion velocities. 
Choosing rise times of 14.0-15.5 d, most weaker lines were 
fitted nicely, while in stronger lines such as Sill A6355 or 
Ca II H&K there was still a lack of absorption at low velocity. 
While weaker lines mostly form near the photosphere, strong 
lines have considerable strength also in layers well above. 
Thus, their centroid typically shows a larger blueshift, and 
their emission, centred at the rest wavelength, is more pro- 
nounced. In order to decrease the flux in their red wings, 
Wph had to be reduced to 5000-6000 kms~^, corresponding 
to rise times of 22-20 d. The weaker features, however, then 
appeared at velocities too low in the synthetic spectrum. 

The difficulty to match all line-velocities simultane- 
ously points to some substantial shortcoming in our treat- 
ment of the SN ejecta. Stratifying the ejecta composition 
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Figure 19. Velocity measured from the minimum of the Sill 
A6355 hne as a function of time. Open green symbols repre- 
sent low-velocity-gradient (LVG) SNe, filled blue symbols high- 
velocity-gradient (HVG) SNe, and red starred s ymbols correspond 
to the FAINT SNe, as defined bv lBenetti et al.i (120051 ). The val- 
ues for HVG and LVG SNe are taken from Benetti et al. and 
IPastorello et al.l l|2007bh . while the data set for the FAINT class 
has been extended and in parts remeasured. 

ijStehle et al.ll2005l ) assuming only little Si and Ca at high 
velocity might cure some of the problems, as it could po- 
tentially remove the blue parts of the respective strong ab- 
sorption lines and reduce their re-emission. Yet, we retained 
our uniform-abundance approach, because no observational 
data earlier than — 6d are available to constrain the addi- 
tional parameters of a stratified-abundance model. Alterna- 
tive explanations for the inconsistencies include a possible 
inadequacy of the W7 density structure, a deviation of the 
total ejected mass from Mch, or 3D effects. 

Nonetheless, based on our model sequence, we con- 
fidently exclude rise times shorter than 14 and longer 
than 22 d. For all spectral models discussed below we as- 
sumed a rise time of 17.0 d. This value is fully compatible 
with the lower limits observationally derived above. At the 
same time, it is slightly shorter than what is favoured for 
Ami5(i?)true = 1.1-SNe, in agreement wit h the trend found 
for SNe la with Ami5(B)truo 1.5 fe.g. iRiess eral]|l999l : 
iKasen fc Wooslevll2007l ). 

7.3 Spectroscopic parameters 

A spectroscopic comparison of SN 2005bl with other under- 
luminous SNe la was performed in Section [5.21 in particular 
in Figs. [51 1111 Here, we concentrate on parameters such as 
7^(Si) (Nug ent et all [l995i ) and V (the average daily rate 
of velocity decrease of Sill A6355 between B-band maxi- 
mum and either the time the Si ii feature dis appears or the 
last a vailable spectru m, whichever is earlier; iBenetti et al.l 
I2OO5I '). With respect to lBenetti et al.l l|2005h . these parame- 



Figure 20. Sin A6355 post-maximum velocity gradient i) 
dBenetti et al.ll2005l) versus Amis (iJ)true- Filled blue squares are 
HVG SNe, open green squares LVG SNe, and filled red stars 
FAINT SNe. 
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Figure 21. Si line-depth ratio 7^(Si)max jNugent et al .Il995h as a 
function of SN decline rate Ami5(_B)truo- The strong correlation 
between these two parameters is evident. Symbols as in Fig. 1201 



ters have been remeasured for s ome FAINT SNe usin g ad- 
ditional data, and SNe 1998de (|Matheson et al.ll2007l ') and 
2005bl have been added. 

Fig. [19] shows the velocity evolution of the Si 11 A6355 
line of SN 2005bl and a large set of comparison objects. 
Once more, SN 2005bl is extremely similar to SNe 1991bg, 
1997cn and 1999by. Before maximum light, the Si velocities 
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Figure 22. Time-evolution of the iNueent et al. ratio 
7^(Si) before and near maximum light. The same symbols as 
in Fig. 1191 have been used. Interpolating lines between adjacent 
points have been drawn to guide the eye. 



of these SNe are comparable to those of the slower members 
of Benetti et al.'s LVG group, but the velocity decrease is 
much faster. Hence, after maximum hght, they have clearly 
the lowest expansion velocities. SN 1998de, photometrically 
nearly a twin of SN 2005bl, has noticeably higher expan- 
sion velocities at early phases. Whether this is caused by 
a high- velocity component (jMazzali et al.|[2005t ). by differ- 
ences in the density structure, or by a higher kinetic en- 
ergy is difflcult to decide without detailed modelling. High- 
velocity features could provide a natural explanation for the 
observed differences without the need to change the explo- 
sion energetics. On the other hand the Sill A6355 line of 
SN 1998de looks symmetric, and the velocities are higher 
than in other 91bg-like SNe not only at the earliest phases, 
but also near maximum hght, when in most other SNe the 
high- velocity feature s have disappeared (Mazzali ot al. 2005; 
iGaravini et al.1 120071 '). Finally, the transitional SNe 1993H 
and 1986G have higher velocities than SN 2005bl and the 
other 91bg-like SNe (except SN 1998de), more similar to the 
LVG group. Thei r comparatively shal low post-maximum ve- 
locity decrease v (jBenetti et al.ll2005f ) also shows their prox- 
imity to LVG SNe. 

In the V vs. Ami5(i3)truc plane, FAINT SNe seem to 
form a separate cluster (Fig. I20p . However, this impression 
may be caused by the lack of SNe with Ami5(i3)truo be- 
tween 1.5 and 1.7 in the Benetti et al. sample. Whether or 
not this gap is real is not easy to decide given the poor 
statistics for rapidly-declining SNe la. With respect to v 
alone, FAINT SNe are very heterogeneous, with values rang- 
ing from those typical of LVG SNe to the highest values 
ever recorded in a SN la (for SN 1998de). The mean i) of 
95 ± 27 (statistical) kms~^d~^ is similar to that reported 
by Benetti et al. for HVG SNe (97± 16 kms'^ d"^). Within 
the FAINT group there may be a tendency for higher v with 



larger Ami5(B)truo, but there are too few objects and their 
range in Amis (S)truo is too small to postulate a correlation. 

Concerning 7?.(Si)max, the value of 7J(Si) at _B-band 
maximum, we confirm the observed trend of larger values 
for faster decliners (Fig. I2ip . The origin of this behaviour 
was the subject of a number of studies, as an explanation 
is not straightforward. iHachinger. Mazzali fc Benettil (I2OO6I ) 
have shown that the change in 7?.(Si)max with Ami5(_B)truo 
is almost solely caused by a variation in the strength of the 
Sill A5972 line, which is somewhat counter-intuitive as the 
excita tion of this line should b e favoured by higher tempera- 
tures. iGarnavich et al.l (|2004 ) suggested that in underlumi- 
nous SNe la the feature at ~ 5800 A might be dominated by 
Ti lines. However, from the synthetic spectra sho wn in Sec- 
tion | 6 .2I we cannot confirm this option. Instead. iHachingerl 
12007? and in prep.) demonstrate that the observed tendency 
is a combined ionisation and excitation effect, natural rather 
than unexpected. 

The pr e-maximum time-evo lution of 7?.(Si), which was 
also used bv lBenetti et al.1 (|2005l ) to distinguish between the 
LVG and HVG groups, cannot be studied equally well for 
the FAINT SNe because of the lack of sufficiently early data 
(more than 5d before maximum). In fact, only before day 
—5 do LVG and HVG SNe show significant differences in 
their behaviour, 7?.(Si) decreasing with time for HVG SNe, 
and being constant or increasing slightly for LVG SNe. Near 
maximum light, the evolution seems to be fairly fiat for 
all SNe la, including the FAINT SNe which merely exhibit 
larger absolute values of 7?.(Si) (see Fig. I22|l . However, es- 
pecially the Sill A5972 line may suffer from variable con- 
tributions of other elements, notably NalD at later phases 
when the spectra become cooler. This limits the conclusive 
power of all studies related to the evolution of 7?.(Si) after 
maximum light, as already at day -|-5 NalD is visible as 
a shoulder in the blue wing of the Sill A5972 feature (cf. 
Section 15. If) . 



7.4 Element abundances from synthetic spectra 

Spectrum synthesis calculations offer an insight into the 
chemical composition of the ejecta, making it possible to 
trace the nuclear reactions that took place during explosive 
burning. However, because of the limitations and possible 
shortcomings of the models mentioned in Section |6] (i.e., no 
abundance stratification, use of a Chandrasekhar-mass pro- 
genitor and of the W7 density distribution), the exact num- 
bers in Table [7] should be taken with caution. Therefore, we 
confine ourselves to analysing clear trends, which should be 
robust with respect to refined models. 



Unburned material, carbon detection 

As for all rapidly-declining SNe la the spectrum of SN 2005bl 
is characterised by strong O I A7774. The strength of this 
line is largely a temperature effect, as the comparatively 
low temperatures encountered in 91bg-like SNe la result in 
a fair fraction of neutral oxygen besides the singly ionised 
state which normally dominates the ejecta. Nevertheless, 
it also signals the presence of a significant amount of un- 
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Figure 23. Identification of carbon in the spectra of SN 2005bl. A 
comparison of the —6 d spectrum of SN 2005bl and the —4 d spec- 
trum of SN 1999by shows a flux deficit at 6400 A in SN 2005bl, 
hkely caused by C II A6580. The insert shows synthetic spectra for 
day —6, one with the composition reported in Table [7] including 
4.5% of C, the other without any C. 

burned material in the SN ejectalHl This is fully consistent 
with the often-anticipated low burning efficiency of these 
objects (from light curves and nebular spectra one can in- 
fer that the amount of synthesised ^''Ni is ^Q.IMq, which 
is at most a quarter o f that in normal- luminosity SNe la, 
IStritzinger et al]|2006l ). Direct information on the oxygen 
abundance is difficult to derive from synthetic spectra, since 
O I A7774 is the only strong feature of this element at op- 
tical wavelengths, and is heavily saturated in all our mod- 
els. Hence, moderate changes in the oxygen content have 
no impact on the strength of this feature. Additionally, our 
abundance estimate should be considered an upper limit, as, 
guided by nuclear reaction network calculations, we ascribe 
the entire mass without observable signatures to oxygen. 

Applying this strategy, we find oxygen mass fractions 
ranging from > 90% at day —6 to ~ 80% at day +5. With a 
W7 density profile and photospheric velocities of 7500 and 
6000 kms~^, about 50 and 60% of the ejecta mass are lo- 
cated above the photospheres at days —6 and +5, respec- 
tively. Hence, if a Chandrasekhar-mass explosion with W7 
density profile provides an acceptable model for 91bg-like 
SNe la, at least the outer 50% of the ejecta of SN 2005bl 
and other members of this group are entirely dominated by 
unburned material. These numbers are significantly differ- 
ent from the results obtain ed by modelling normal SNe la 
(cf., e.g.,|stehieelaD[200i). They are also in disagreement 
with the picture sketched in the Zorro plot l|Mazzali et al.l 
I2OO7I '). i.e. that all LVG and FAINT SNe la have about the 
same amount of unburned material, and that the observed 
spectrophotometric sequence among SNe la has its origin 
only in a variable ratio of ^^Ni to IMEs. 

Another indication of the low burning efficiency in 
SN 2005bl is the likely detection of unprocessed carbon in its 

^■^ In this section we refer to C and O as unburned, although 
— depending on the initial C-to-O ratio of the white dwarf - a 
significant fraction of the O may have been produced during the 
explosion by incomplete C burning. 



early spectra. Although the S/N of the —6 and — 5d spec- 
tra is low, a clear flux deficit redwards of the Sill A6355 
absorption can be discerned compared to coeval spectra of 
SN 1999by (Fig.[23l). 

Synthetic spectra enable us to investigate the effect of 
including different chemical species on this region. There are 
only few lines which could potentially contribute to the ob- 
served absorption 13 Including hydrogen did not improve the 
model, since the resulting Ha absorption was at too short 
a wavelength, and Ha emission deteriorated the fit around 
6500 A. A clearly better match to the observed spectrum was 
obtained introducing carbon, as C 11 lines not only cropped 
the flux peak at 6400 A through Cll A6580 (Fig. [Ml insert), 
but also improved the fit in other regions, most notably 
around 7000 A (through Cll A7234). Therefore, we consider 
the presence of a few percent of carbon in the ejecta at days 
—6 and —5 very likely. 

The temperature and density conditions for carbon 
burning are more relaxed than those for oxygen burning, 
so that carbon is burned much more completely to heav- 
ier elements than oxygen. In fact, only in some - usually 
very early (~— lOd) - optical spectra of SNe la have clear 
signatures o f carbon ever been found ( most prominently in 
SN 2006gz, iHicken et all l2007l : cf. also [Branch et al] |2007| 
for an overview). None of these SNe was 91bg-like, although 
the lower burning efficiency in these objects might favour 
the presence of unprocessed carbon. The lack of spectra 
taken earlier than 5d before maximum light may have so 
far prevented the detection of carbon. In SN 2005bl, how- 
ever, we are confident that signatures of carbon are present 
less than one week before maximum. This constitutes one 
of the few detections of C 11 in SN la spectra at such a rel- 
atively late epoch; previou s cases include SNe 1996X and 
2006D ([Thomas et al.ll2007l ). 

Abundance of Fe-group elements 



The Fe-group abundance of SN 2005bl of ~ 0.01% at day 
—6 and ~ 0.4% at day -f 5 is about two orders of mag- 
nitude lower than in normal-luminosity SNe la at com- 



parable epochs (see, e.g., 


Stehle et all 20051: Kotak et al.[ 


[2005l:[Ehas-Rosa et al.|[200e 


: Altavilla et al. 2007). At a first 



glance this looks perfectly consistent with the low luminos- 
ity and hence the low production of ^^Ni and other NSE- 
elements in a 91bg-like SN la. 



However, closer examination reveals that the Fe con- 
tent of the pre-maximum models is significantly sub-solar. 
At — 6d the best-fitting model contains no Ni or Co, and 
only 0.010 % of stable Fe. Acceptable results were obtained 
with Fe abundances up to 0.025%, but also without any 
Fe (see Fig. I24|l . Solar Fe abundances (the F e mass fraction 
of the Sun's photosphere is about 0.115 %. [Asplund et al.[ 
(2005i ) can safely be ruled out, as with this amount of Fe the 



Conversely, because of the relative lack of strong lines, the re- 
gion around 6400 A proved to be a preferential window for flux 
redistributed from shorter wavelengths to escape. All elements 
which create opacity in the blue and UV part of the spectrum 
(Fe-group elements, Ti, Cr, other metals) rather led to an in- 
crease of the 6400 A flux excess in the model. This in turn yielded 
constraints on the metal abundance. 
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Figure 24. Sequence of synthetic spectra, varying the mass frac- 
tion of Fe between to 0.1% (the latter roughly represents so- 
lar Fe abundance). The observed but smoothed — 6d spectrum is 
shown in black. The parameters of the best-fitting spectrum (with 
X(Fe)= 0.01 %) are summarised in Table [7] The other spectra 
have been obtained by changing only the amount of stable Fe at 
the expense of O, the most abundant element, without attempting 
to optimise the fit by adapting other parameters. Hence, better 
fits might be feasible also for the models with X(Fe)^ 0.01 %. 



quality of the fit deteriorates a lot around 5000 A. Now, the 
Fe content deduced not only accounts for ^''Pe synthesised 
during the explosive burning (the ^^Ni decay chain has too 
long decay times to give a relevant contribution at these 
epochs), but also for Fe already present in the progenitor 
star before the explosion. Hence, a sub-solar Fe content of 
the precursor star might be inferred. On the other hand, 
SN 2005bl exploded in an elliptical galaxy, a presumably 
not metal-poor environment. If the chemical composition of 
the precursor did not deviate significantly from the average 
material in NGC 4070, it rather should have had a fairly 
high Fe abundance, and nuclear burning should have fur- 
ther enriched the ejecta with Fe-group elements, although 
this may be confined to the inner shells. 

There may be two ways to explain the small amount of 
Fe required to fit the early-time spectra of SN 2005bl. One is 
to assume that Fe was not distributed homogeneously in the 
progenitor star, but concentrated towards the centre. A cer- 
tain degree of gravitational settling of heavier elements can 
indeed be expected. Estimating whether or not this effect 
is sufficient to yield Fe abundances in accordance to those 
inferred from the synthetic spectra would require a model of 
gravitational diffusion inside the white dwarf, and accurate 
constraints on the mass accretion rate from the companion 
and the composition of the accreted material. Note, however, 
that assuming Wph = 7500 kms~^ and a W7 density profile 
the photosphere at day —6 is not located in the outermost 
layers, but half the way down the ejecta in mass coordinates. 

Alternatively, the progenitor of SN 2005bl might have 
had a metallicity truly lower than the average metallicity of 
an elliptical galaxy. This could be explained if the progen- 
itor was rather old (~ 10 Gyr), formed before the ISM in 
NGC 4070 was metal-enriched through recurrent cycles of 
stellar birth, mass loss and death. If this scenario held true, 
it would provide one of the first direct constraints on the 



lifetime of an individual SN la progenitor. Given the spec- 
troscopic similarity of SN 2005bl to SNe 1991bg, 1997cn, 
1998de and 1999by, and given that most of the latter SNe 
also exploded in early-type, supposedly metal-rich galaxies 
(see Table |8]), this could point towards a very long-lived pro- 
genitor population for all underluminous, 91bg-like SNe la. 

A speculation could be that 91bg-like SNe la descend 
from white-dwarf mergers (DD scenario). White-dwarf bi- 
nary systems should meet the requirement of long life- 
times because of the low efficiency of angular-momentum 
loss through gravitational-wave emission. In fact, in a study 
of the delay-time distribution of SNe la, iGreggiol (|2005l ) 
identified initially wide white-dwarf pairs as the only pro- 
genitor model whose delay-time distributi on does not drop 
drama tically beyond ~ 8 Gyr. Models of iKobavashi et al.1 
(|l998i) lend further support to the DD scenario, as they pre- 
dict that no SNe la should occur in SD systems for [Fe/H] 
^ —1, a limit which may be exceeded by SN 2005bl. Fi- 
nally, from s p ectro polarimetric observations of SN 1999by 
iHowell et al.l (|200ll ) inferred a degree of asphericity much 
higher than in ordinary SNe la. On this basis they identified 
rapidly rotating white dwarfs and merging white-dwarf bina- 
ries as the most promising progenitor candidates for 91bg- 
like SNe la. A different progenitor system with respect to 
ordinary SNe la (for which the SD scenario consisting of a 
white dwarf and a non-degenerate companion is favoured), 
and in particular the possibility that the mass is different 
from Mch , could account for observed differences in the pho- 
tometric and spectroscopic properties of 91bg-like SNe. 



8 CONCLUSIONS 

We have presented and analysed optical photometric and 
spectroscopic data of the underluminous Type la SN 2005bl 
from one week before to two months after maximum light 
in B. In this entire interval the evolution of SN 2005bl is 
substantially different from that of slow or intermediate de- 
cliners, but very similar to that of the prototypical fast de- 
cliners SNe 1991bg and 1999by. At peak, SN 2005bl appears 
slightly more luminous than the latter, but the difference is 
marginal and within the error bars, which for SN 2005bl are 
dominated by uncertainties in the amount of host-galaxy ex- 
tinction. With respect to quantities like the Ami5(_B)truc of 
1.93, the ejecta velocities inferred from the minimum of the 
Sill A6355 line, and temperature indicators such as 7?.(Si) 
or the time-evolution of various colour indices, SN 2005bl is 
in fact - within the measurement uncertainties - a clone of 
SNe 1991bg and 1999by. 

We confirm that Ami5(B)truo is a good luminosity in- 
dicator also for underluminous SNe la, with a steeper de- 
pendence of the peak magnitudes on Ami5(_B)truo than 
for slow and intermediate decliners. From linear fits we 
obtain , = 6.83 ± 0.32 and 4.33 ± 0.31 for 

d(Ami5(S)truo) 

the B and V band, respectively, considering objects with 

Ami5(B)true ^ 1.69. 

We conducted an analysis analogous to that of 
iBenetti et al.l (j2005l ). focussing on the behaviour of the 
FAINT subclass and trying to identify trends within this 
group. In addition to the SNe in Benetti et al. we not only 
added the measurements of SN 2005bl, but also inc luded 
SN 1998de (jModiaz et al.ll200ll : [Matheson et al.ll2007l ) and 
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updated some numbers on the basis of better data avail- 
ability. We find a correlation between Ami5(_B)truc and 
7^(Si)max, nicely extending the observed linear trend among 
HVG and LVG SNe to larger decline rates. Because of the 
lack of early data, the time-evolution of 7?.(Si) at phases 
where differences between the LVG and IfVG groups be- 
come evident cannot be investigated. The Sill line velocities 
of FAINT SNe are mostly similar to those of LVG SNe near 
maximum light, but significantly lower a few weeks later. 
There is a large dispersion in post-maximum velocity gra- 
dients i} within the FAINT class, with numbers comparable 
to those of HVG SNe, but exceeding their range in both 
directions. 

Synthetic spectra for SN 2005bl at four epochs before 
and soon after maximum light have been computed using a 
ID Monte Carlo code. Despite uncertainties and shortcom- 
ings such as the poorly-constrained rise time, the possible in- 
adequacy of the W7 density profile for 91bg-like SNe la, and 
the fact that computations were performed without chem- 
ical stratification, a number of interesting results were ob- 
tained. The presence of carbon in the —6 and — 5d spectra, 
indicated by a visual comparison with SN 1999by at similar 
epoch, was confirmed, and an overall low burning efficiency 
was established. NSE elements, but also IMEs, are signifi- 
cantly less abundant than in ordinary SNe la; instead, most 
of the ejecta above 6000 kms~^ are made of oxygen, and 
are hence either unburned or the result of highly incomplete 
carbon burning. Furthermore, in the —6 and — 5d spectra 
only traces of Fe are found, one order of magnitude less than 
expected for unprocessed material with solar composition. 
Possible implications of this low metallicity on the nature of 
the progenitor were discussed in Section [7.41 

It is presently very difficult to decide whether 91bg-like 
SNe la form a separate group, distinct from ordinary SNe la 
and possibly descending from different progenitors or ex- 
plosion mechanisms, or whether they are just the extreme 
end of a continuous distribution of objects, with their char- 
acteristic appearance owing to lower temperatures in con- 
junction with a different chemical composition, both caused 
by the overall low burning efficiency. To tackle this prob- 
le m will require a new statistical analysis similar to that 
of iBenetti et al.l (|2005h once a bigger sample of rapidly- 
declining SNe la will be available, and more detailed spectral 
modelling of the data sets already published. 
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